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Abstract 

Interference  and  jamming  signals  are  a  serious  concern  in  an  operational  military  communica¬ 
tion  environment.  This  thesis  examines  the  utility  and  performance  of  combining  adaptive  temporal 
filtering  with  adaptive  spatial  filtering  (i.e.  adaptive  beamforming)  to  improve  the  signal-to-jammer 
ratio  (SJR)  in  the  presence  of  narrowband  and  wideband  interference.  Adaptive  temporal  filters  are 
used  for  narrowband  interference  suppression  while  adaptive  beamforming  is  used  to  suppress  wide¬ 
band  interference  signals.  A  procedure  is  presented  for  the  design  and  implementation  of  a  linear 
constraints  minimum  variance  generalized  sidelobe  canceler  (LCMV-GSC)  beamformer.  The  adap¬ 
tive  beamformer  processes  the  desired  signal  with  unity  gain  while  simultaneously  and  adaptively 
minimizing  the  output  due  to  any  undesired  signal.  Using  the  LCMV-GSC  beamformer  with  a  least 
mean  squares  (LMS)  adaptive  algorithm,  it  was  shown  that  the  tapped  delay  line  (TDL)  adaptive 
antenna  array  is  more  effective  for  the  suppression  of  wideband  jammer  suppression  than  the  linear 
array  sensors  (LAS)  adaptive  antenna  array.  Also  a  new  technique  for  adaptive  beamforming  is  pre¬ 
sented  which  improves  wideband  interference  suppression  in  a  frequency-hopped  environment.  The 
output  SJR  improvement  for  the  new  technique  compared  to  the  conventional  technique  is  as  much 
as  15dB.  Sometimes,  multipath  signals  and  jammers  generated  by  a  smart  enemy  are  correlated 
with  the  desired  signal  which  destroys  the  traditional  beamformer’s  performance.  After  perform¬ 
ing  a  spatial  smoothing  technique,  adaptive  beamforming  can  also  be  effective  in  suppressing  the 
jamming  signals  that  are  highly  correlated  with  the  desired  signal. 
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Interference  Suppression  for  Spread  Spectrum  Signals 
Using  Adaptive  Beamforming  and  Adaptive  Temporal  Filter 

1.  Introduction 

1.1  Background 

The  final  goal  of  communication  is  to  transmit  information  without  error  from  a  transmit¬ 
ter  to  a  receiver.  In  the  real  world,  there  are  several  interference  sources  with  various  powers 
and  bandwidths  to  prevent  this  goal.  These  interfering  signals  might  arise  from  intentional  jam¬ 
ming,  multipath,  multiple  users  on  the  same  bandwidth,  or  a  variety  of  other  sources.  In  a  noisy, 
crowded,  or  hostile  environment,  a  number  of  techniques  may  be  implemented  in  order  to  increase 
the  anti-jam  characteristics  of  a  communication  system.  Some  methods  include  spread  spectrum 
modulation,  interference  suppression  and  spatial  discrimination  techniques. 

The  objective  of  this  research  is  to  examine  the  effectiveness  of  adaptive  algorithms  for  in¬ 
terference  suppression  in  spread  spectrum  communication.  It  is  desirable  to  achieve  temporal  and 
spatial  filtering  using  adaptive  algorithms  to  maintain  an  acceptable  signal-to-noise  ratio  (and 
hence,  probability  of  bit  error)  in  any  kind  of  interference  such  as  intentional  jamming,  multipath, 
and  co-channel  interference. 

1.1.1  Spread  Spectrum  Modulation.  Spread  Spectrum  techniques  have  an  inherent  ability 
to  suppress  interference  (Peterson  at  el,  1995).  These  are  modulation  schemes  that  produce  a 
spectrum  with  a  much  wider  bandwidth  than  that  of  the  information  bearing  signal.  By  increasing 
the  bandwidth  over  which  signal  energy  is  contained,  the  energy  density  of  the  signal  is  reduced 
and  hidden  in  the  noise  levels  so  that  unintended  receivers  cannot  detect  it.  Hence  the  signals  have 
a  low  probability  of  intercept.  Upon  reception  of  the  desired  signal  and  undesired  jammer  signal 
by  the  Spread  Spectrum  (SS)  receiver,  the  coding  process  is  reversed  and  the  signal  energy  density 
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is  returned  to  its  original  level  while  the  jammer  energy  density  is  reduced.  Hence  the  receiver  is 
resistant  to  interference  and  jamming. 

Two  common  spread  spectrum  schemes  are  direct  sequence  (DS)  and  frequency  hopping 
(FH).  Direct  Sequence  is  achieved  by  combining  a  binary  message  sequence  and  a  higher  rate 
pseudorandom  binary  sequence.  A  pseudorandom  binary  sequence  code  has  a  much  higher  rate 
than  the  original  digital  data  which  expands  the  bandwidth  beyond  that  of  the  original  information 
bandwidth.  Frequency  hopping  is  achieved  by  changing  the  frequency  of  the  carrier  periodically. 
Typically,  each  carrier  frequency  is  chosen  from  a  set  of  2^  frequencies  which  are  spaced  apart 
approximately  the  width  of  the  data  modulation  bandwidth.  The  frequency  is  translated  to  one 
of  2^  frequency  hop  bands  by  the  FH  modulator.  The  processing  gain  is  2^.  The  spreading  code 
in  this  case  does  not  directly  modulate  the  carrier  but  is  instead  used  to  control  the  sequence  of 
carrier  frequencies. 

1.1.2  Interference  Suppression.  The  inherent  processing  gain  of  a  spread  spectrum  signal 
will,  in  many  case,  provide  the  system  with  a  sufficient  degree  of  interference  rejection  capability. 
However,  if  the  combined  interference  signal  power  relative  to  the  desired  signal  power  exceeds  the 
spread  spectrum  processing  gain,  additional  filtering  is  required.  If  the  interference  is  relatively 
narrowband  compared  with  the  bandwidth  of  the  spread  spectrum,  then  interference  suppression 
by  the  use  of  notch  filters  often  results  in  an  improvement  in  system  performance. 

There  are  two  techniques  for  building  notch  filters.  The  first  technique  uses  a  transversal  filter 
in  the  time  domain  (Milstain  and  Das,  1980;  Saulnier,  1990).  The  system  can  be  made  adaptive 
by  using  a  tapped  delay  line  with  variable  tap  weights.  These  tap  weights  can  be  adapted,  for 
example,  by  using  the  well-known  least-mean-square  (LMS)  algorithm. 

The  second  technique  is  that  of  transform  domain  processing.  A  notch  filter  is  implemented 
by  Fourier  transforming  the  received  waveform,  applying  some  type  of  signal  processing  algorithm, 
and  then  inverse  Fourier  transforming  the  signal  back  to  the  time  domain.  One  type  of  signal 
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processing  algorithm  excises  the  frequency  bin  which  exceeds  a  specific  threshold  (Milstain  and 
Das,  1980).  Another  way  is  to  use  adaptive  transform  domain  filtering  (Saulnier,  1992).  This 
type  of  suppressor  combines  features  of  the  time  domain  adaptive  filter  and  the  transform  domain 
excisor.  The  transform  domain  technique  works  more  effectively  than  the  time  domain  technique 
for  all  interference  bandwidth.  Transform  domain  techniques  are  also  more  robust  since  they  can 
handle  multiple  narrowband  jammers  simultaneously. 

1.1.3  Spatial  Discrimination  Technique.  Spatial  discrimination  techniques  include  adap¬ 
tive  null  steering  antennas  and  high-gain  directional  antennas  (Widrow  et  al,  1982).  Together, 
they  place  nulls  of  the  receiver’s  antenna  pattern  in  the  direction  of  enemy  jammers  to  avoid  front 
end  saturation  of  the  receiver  and  to  increase  the  gain  in  the  direction  of  the  desired  signal.  In  an 
adaptive  array,  the  phase  and  amplitude  of  the  signal  at  each  receiving  antenna  element  is  weighted 
and  the  resulting  signals  summed  to  produce  the  array  output.  The  values  of  the  element  weights 
are  determined  by  an  algorithm  that  can  act  to  steer  nulls  in  the  direction  of  interfering  signals 
(Fig.  1.1). 

1.2  Problem  Statement  and  Objective 

Any  kind  of  interference  signals  are  a  serious  concern  in  an  operational  military  communica¬ 
tion  environment.  Many  researchers  have  developed  lots  of  techniques  for  narrowband  interference 
suppression,  i.e.  adaptive  filter  in  time  and  transform  domain  and  excisor.  In  the  case  of  wideband 
interferences,  those  techniques  can’t  work,  because  the  desired  signal  is  also  suppressed  when  they 
suppress  the  wideband  interferences.  Until  now,  the  adaptive  beamforming  is  known  as  the  way  to 
suppress  the  wideband  interferences.  However,  the  wideband  interferences  suppression  remains  an 
important  research  topic.  Furthermore,  multipath  signals  are  correlated  to  the  desired  signal  and 
the  smart  jammer  uses  the  correlated  signal  with  the  desired  signal  to  prevent  our  communication. 
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The  conventional  adaptive  beamforming  is  totally  destroyed  in  the  presence  of  correlated  signals. 
The  correlated  interference  is  a  serious  problem  in  the  traditional  adaptive  beamformer. 

This  research  combines  a  notch  filter  technique  and  spatial  discrimination  technique  which 
are  two  of  the  most  powerful  methods  for  providing  interference  signal  suppression.  In  other  words, 
by  achieving  adaptive  temporal  and  spatial  filtering  using  adaptive  algorithm,  the  system  maintains 
the  reasonable  signal  power  to  jammer  power  ratio  (SJR). 

1.3  Assumption 

The  research  assumes  that  the  received  signal  is  composed  of  a  sum  of  jamming  and  desired 
signals.  It  is  also  assumed  that  code  synchronization  is  maintained  at  the  communication  receiver, 
and  thus,  there  is  no  propagation  delay.  Finally,  unless  explicitly  mentioned,  the  desired  signal  and 
jamming  signal  are  uncorrelated  each  other.  The  antenna  elements  are  uniformly  distributed  and 
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the  distance  between  elements  is  a  half  wavelength  of  the  desired  signal  center  frequency  in  order  to 
prevent  the  grating  lobe.  Those  assumptions  are  reasonable  for  real  situation  and  hardware,  thus 
they  doesn’t  affect  the  result  of  simulation. 

1.4  Scope 

The  desired  signal  considered  in  this  research  has  wide  bandwidth.  Specifically,  frequency 
hopped  signals  are  considered  because  FH  (Frequency  Hopping)  signals  and  antenna  arrays  are 
closely  related.  Each  antenna  element  is  omni-directional.  Jammer  models  consist  of  both  narrow 
and  wide  bandwidth.  The  adaptive  algorithm  used  in  both  schemes  is  the  LMS  (Least  Mean 
Squares)  algorithm.  This  algorithm  is  the  best  known  and  most  easily  implemented  algorithm 
which  implement  an  iterative  solution  to  the  Wiener-Hoff  equation  without  making  use  of  any  a 
priori  statistical  information  about  the  received  signal.  The  Wiener-Hoff  equation  determines  the 
optimal  tap  weight  settings  for  the  transversal  filter. 

The  class  of  adaptive  beamformer  discussed  is  LCMV-LMS  (Linearly  Constrained  Minimum 
Variance  -  Least  Mean  Squares)  beamforming  which  constrains  the  gain  and  phase  in  the  desired 
signal  direction  and  minimizes  the  output  due  to  the  undesired  signals  (Griffiths  and  Jim,  1982). 
For  LCMV  beamforming,  a  GSC  (Generalized  Side-lobe  Canceller)  structure  is  used  because  this 
structure  is  useful  for  implementation  and  analysis  of  LCMV  beamforming. 

1.5  Material  and  Equipment 

All  system  models  and  simulations  were  developed  using  the  version  4.1  Matlab  simulation 
software  developed  by  The  Math  Works,  Natric,  Massachusetts.  The  software  was  run  on  a  Sun 
Workstation. 
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1.6  Thesis  Organization 


Chapter  II  of  this  thesis  highlights  a  literature  review  of  adaptive  algorithm,  adaptive  beam¬ 
forming,  and  adaptive  signal  processing  in  spread  spectrum  modulation.  Chapter  III  presents  a 
basic  introduction  into  the  concept  of  adaptive  signal  processing.  It  compares  the  temporal  fil¬ 
ter  versus  spatial  filter.  It  then  discusses  the  theory  of  adaptive  temporal  filters  and  adaptive 
beamforming.  Moreover,  the  adaptive  beamforming  in  the  presence  of  correlated  signals  is  pre¬ 
sented.  Chapter  IV  provides  the  simulation  results  related  to  jamming  suppression  and  especially 
emphasizes  the  jamming  suppression  in  frequency-hopped  environment.  In  addition,  it  provides 
the  performance  of  a  spatial  smoothing  technique  for  suppression  of  correlated  jammers.  Chapter 
V  summarizes  the  thesis,  states  conclusions  based  on  the  simulation  results,  and  provides  recom¬ 
mendations  for  future  research. 
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II.  Literature  Review 


This  chapter  reviews  the  published  literature  on  the  interference  rejection  problems  as  it 
relates  to  adaptive  signal  processing.  Adaptive  filtering  algorithms  and  adaptive  beamforming  for 
wideband  interferences  and  for  the  specific  case  spread  spectrum  signals  are  reviewed. 

2.1  Adaptive  Filtering  Algorithm 

The  first  studies  of  minimum  mean-square  estimation  in  stochastic  processes  may  be  traced 
back  to  the  late  1930s  and  early  1940s.  Wiener  formulated  the  continuous  time  linear  prediction 
problem  and  derived  an  explicit  formula  for  the  optimum  predictor.  He  also  considered  the  filtering 
problem  of  estimating  a  process  corrupted  by  an  additive  noise  process.  The  explicit  formula  for  the 
optimum  estimate  required  the  solution  of  an  integral  equation  known  as  the  Wiener-Hoff  equation 
(Wiener  and  Hoff,  1931) 

In  1947,  Levinson  formulated  the  Wiener  filtering  problem  in  discrete  time.  The  Wiener-Hoff 
equation  can  be  expressed  in  matrix  form  as  follows, 

Rwo  =  p  (2.1) 

where  Wq  is  the  tap  weight  vector  of  the  optimum  Wiener  filter  in  the  form  of  a  transversal  filter, 
which  is  shown  in  Fig.  2.1.  R  is  the  correlation  matrix  of  the  tap  inputs  u(7i),  and  p  is  the 
cross-correlation  matrix  of  the  tap  inputs  u(7r)  and  the  desired  response  d{n)  shown  in  Fig.  2.1. 
(Levinson,  1947). 

The  simple  adaptive  filtering  least-mean-square  (LMS)  algorithm  emerged  as  a  algorithm  for 
the  operation  of  adaptive  transversal  filters  in  the  late  1950s.  The  LMS  algorithm  was  devised  by 
Widrow  and  Hoff  in  1959  in  their  study  of  a  pattern  recognition  scheme  known  as  the  adaptive 
linear  threshold  logic  element.  The  LMS  algorithm  is  a  stochastic  gradient  algorithm  in  that  it 
iterates  each  tap  weight  in  the  transversal  filter  in  the  direction  of  the  gradient  of  the  squared 
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Figure  2.1  Transversal  Filter 

amplitude  of  an  error  signal  e(n)  with  respect  to  that  tap  weight.  (Widrow,  1970).  Fig.  2.2  shows 
a  block  diagram  of  the  LMS  method.  The  weights  are  adjusted  or  adapted  in  response  to  the  data 
received  at  the  transversal  filter.  Changing  the  frequency  response  is  as  simple  as  changing  the 
weights.  The  weights  in  an  adaptive  beamformer  are  controlled  by  the  LMS  algorithm. 

Godarn  made  another  major  contribution  using  Kalman  filter  theory  to  propose  a  new  class 
of  adaptive  filtering  algorithms  for  obtaining  rapid  convergence  of  the  tap  weights  of  a  transversal 
filter  to  their  optimum  settings  (Godard,  1974).  This  algorithm  is  referred  to  in  the  literature 
as  the  Kalman  algorithm  or  Godarn  algorithm.  The  Kalman  algorithm  is  closely  related  to  the 
recursive  least-squares  (RLS)  algorithm  that  follows  from  the  method  of  least  squares  (Plackett, 
1950).  The  Kalman  or  RLS  algorithm  usually  provides  a  much  faster  rate  of  convergence  than 
the  LMS  algorithm  at  the  expense  of  increased  computational  complexity.  The  desire  to  reduce 
computational  complexity  to  a  level  comparable  to  that  of  the  simple  LMS  algorithm  prompted 
the  search  for  a  computational  efficient  RLS  algorithm. 

2.2  Adaptive  Beamforming 

2.2.1  General  Aspect.  The  main  object  of  adaptive  beamforming  is  to  put  a  null  in  the 
jammer  direction  using  an  adaptive  filtering  algorithm.  The  technology  of  adaptive  beamforming 
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has  been  developed  with  adaptive  filtering  algorithms.  The  initial  work  of  adaptive  beamforming 
may  be  traced  back  to  the  invention  of  the  intermediate  frequency  (IF)  sidelohe  canceler  by  Howells 
in  the  late  1950.  In  his  historical  report,  Howells  described  a  sidelobe  canceler  capable  of  automat¬ 
ically  nulling  out  the  effect  of  a  jammer.  The  sidelobe  antenna  uses  a  primary  (high  gain)  antenna 
and  a  reference  omni-directional  (low  gain)  antenna  to  form  a  two-element  array  with  one  degree 
of  freedom.  (Howells,  1976) 

The  poor  performance  of  the  delay-and-sum  beamformer  (the  transversal  filter)  is  due  to 
the  fact  that  its  response  along  a  direction  of  interest  depends  not  only  on  the  power  of  the  in¬ 
coming  target  but  also  undesirable  contributions  received  from  other  sources  of  interference.  To 
overcome  this  limitation.  Capon  proposed  a  new  beamformer  in  which  the  weight  vector  w(n)  is 
chosen  to  minimize  the  variance  (i.e.,  average  power)  of  the  beamformer  output,  subject  to  the 
constraint  w^{n)s(6)  =  1  for  all  n,  where  s(^)  is  a  prescribed  steering  vector  (Capon,  1969).  This 
constrained  minimization  yields  an  adaptive  beamformer  with  a  minimum  variance  distortionless 
response  (MVDR)  which  processes  the  desired  signal  from  certain  directions  with  specified  gain  and 
phases  using  the  constraints  and  then  minimize  the  output  power  to  suppress  undesired  signals. 
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Figure  2.3  Griffiths-Jim  adaptive  beamformer  structure 


Another  structure,  suggested  by  Frost,  was  an  adaptive  algorithm  designed  to  make  the 
beamformer  cancel  everything  that  does  not  come  from  the  desired  signal  direction  (Frost,  1972). 
A  linear  hard  constraint  prevents  cancelation  from  the  look  direction,  thus  avoiding  the  trivial 
zero  solution.  A  modification  of  LMS  algorithm  to  a  two-step  algorithm  was  used  to  update  the 
coefficients.  The  Frost  adaptive  beamformer  was  a  source  of  inspiration  to  further  development. 

The  Widrow-McCool  beamformer  rearranged  the  Frost  beamformer  as  a  noise  canceler  with  a 
slightly  different  constraint.  The  Griffiths-Jim  adaptive  beamformer  is  a  modified  Widrow-McCool 
beamformer  where  the  constraint  is  included  in  the  structure  (Griffiths  and  Jim,  1982).  In  its  most 
general  form,  the  Griffiths-Jim  structure  replaces  the  constraints  by  a  signal  blocking  matrix  unit 
(Fig.  2.3). 

2.2.2  Adaptive  Beamforming  for  Wideband  Interference  Suppression.  In  1967,  Widrow  et 
al.  suggested  the  use  of  tapped  delay  lines  in  an  adaptive  array  for  wideband  interference  rejection 
(Widrow,  1967).  Figure  2.4  shows  the  common  broadband  beamformer  with  the  tapped  delay  lines. 
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Figure  2.4  TDL  Antenna  Array  for  Wideband  Jammer  Suppression 


This  antenna  is  referred  to  TDL  (Tapped  Delay  line)  array  antenna,  while  the  antenna  shown  in 
Fig.  1.1  is  referred  to  LMS  array  antenna.  The  study  of  tapped  delay  lines  has  since  been  done 

by  several  others.  In  one  study,  Rodgers  and  Compton  compared  the  performance  of  a  two-element 
array  with  two,  three  and  five-tap  delay  lines  using  real  weights  to  that  of  an  array  with  a  single 
complex  weight  behind  each  element  (Rodgers,  1979).  They  assumed  the  space  between  taps  was 
one  quarter  wavelength.  In  1981,  Mayhan,  Simmons,,  and  Cummings  presented  a  mathematical 
analysis  of  how  the  number  of  elements  and  the  number  of  delay-line  taps  affect  the  interference 
cancelation  ratio  as  a  function  of  bandwidth  (Mayhan,  1981).  In  1983,  White  studied  the  trade¬ 
off  between  the  number  of  interfering  signals  and  the  required  number  of  auxiliary  elements  and 
delay-line  taps  in  an  Applebaum  array  (White,  1983).  In  other  words,  he  addressed  the  bandwidth 
capability  of  the  canceler  by  using  tapped  delay  lines  or  simply  by  using  more  auxiliary  elements. 

In  spite  of  these  contributions,  there  was  still  no  simple  explanation  in  the  literature  for 
how  the  number  of  taps  and  the  amount  of  delay  between  taps  should  be  chosen  in  an  adaptive 
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array  to  achieve  a  given  bandwidth  performance.  In  1988,  Compton  addressed  this  question  by 
examining  a  two-element  adaptive  array  with  a  tapped  delay-line  behind  each  element  (Compton, 
1988).  In  the  same  year,  he  also  considered  the  use  of  fast  Fourier  transform  (FFT)  processing 
behind  the  elements  in  adaptive  arrays  (Compton,  1988).  However,  he  concluded  the  performance 
of  FFT  processing  wasn’t  an  improvement  over  tapped  delay  lines.  He  used  the  LMS  algorithm  for 
changing  the  weights  in  both  cases. 

In  1986,  Buckley  and  Griffiths  presented  an  adaptive  broadband  beamforming  structure  which 
added  the  Frost  beamformer  and  TDL  antenna  array.  This  beamformer  employed  a  gradient-based 
weight  adjustment  to  minimize  output  variance  subject  to  a  set  of  J  linear  constraints  on  broadband 
directional  derivatives  in  the  desired  look  direction.  Generalized  sidelobe-cancelling  structure  was 
employed  in  which  a  nonadaptive  beamformer  operates  in  parallel  with  an  adaptive  beamformer 
(Buckley  and  Griffiths,  1986). 

2.3  Adaptive  Signal  Processing  in  Spread  Spectrum 

2.3.1  Direct  Sequence  Spread  Spectrum  with  Adaptive  Filter.  Specific  implementation  of 
an  adaptive  algorithm  for  interference  suppression  in  Direct  Sequence  Spread  Spectrum(DS-SS) 
systems  has  been  investigated  since  the  mid-  to  late-1970’s.  Hsu  and  Giordano  have  laid  down  the 
foundation  in  interference  suppression  in  DS-SS  systems  (Hsu  and  Giordano,  1978).  Their  digital 
whitening  is  accomplished  by  using  a  transversal  filter  whose  coefficients  are  selected  by  either  a 
Wiener  algorithm  or  a  maximum  entropy  algorithm.  Filters  obtained  by  use  of  these  algorithms 
are  evaluated  for  various  jamming  and  signaling  conditions  and  are  found  to  exhibit  comparable 
performance  over  a  wide  range  of  input  signal  and  noise  ratios.  The  Wiener  filter  was  implemented 
recursively  using  a  least-mean-square  criterion  and  a  form  of  Levinson’s  algorithm  in  the  actual 
computation. 

Ketchem  and  Proakis  further  improved  the  results  of  Hsu  and  Giordano  by  combining  the 
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interference  rejection  filter  with  its  matched  filter,  resulting  in  an  overall  linear  filter  having  a  linear 
phase  characteristic.  They  further  defined  the  performance  of  the  SS  receiver  as  measured  in  terms 
of  the  probability  of  error,  which  was  obtained  by  applying  a  Gaussian  assumption  on  the  total 
residual  noise  and  interference  followed  by  Monte  Carlo  simulation.  Ketchum  and  Proakis  also 
investigated  the  size  of  the  interference  suppression  filter  (in  terms  of  number  of  taps)  required  to 
handle  multiple-band  interference.  The  frequency  response  improves  with  increasing  filter  order 
(the  number  of  taps)  in  terms  of  providing  a  deeper  notch  at  the  interfering  tone  frequencies  and 
less  attenuation  in  the  frequency  range  between  notches  (Ketchum  and  Proakis,  1932). 

Saulnier  has  performed  hardware  implementations  of  interference  suppression  filters  using 
three  different  transversal  filter  structures:  charge  transfer  device  (CTD),  digital  filter  techniques, 
and  surface  acoustic  wave  (SAW)  device  (Saulnier  et  al,  1984;  Saulnier  et  al,  1985;  Saulnier,  1990). 
These  three  implementations  were  of  the  estimation-type  adaptive  filter.  A  CTD-based  adaptive 
filter  was  designed  using  the  Widrow-Hoff  LMS  algorithm  to  implement  the  adaptive  filter  archi¬ 
tecture  consisting  of  only  two  multipliers  (Saulnier  et  al,  1984).  This  hardware  simplification  was 
achieved  through  the  use  of  a  burst  processing  technique.  Saulnier,  Das,  and  Milstein  followed  this 
work  with  an  identical  setup  with  a  digital  hardware  implementation  instead  of  the  CTD  (Saulnier 
et  al,  1985).  They  recently  performed  a  hardware  experimentation  using  a  SAW-based  adaptive 
filter  to  perform  the  interference  suppression  (Saulnier,  1990).  The  primary  advantage  of  using  the 
SAW-based  transversal  filter  is  the  availability  of  higher  bandwidths  and  operating  frequencies  (i.e. 
correlation  of  the  DS  signal  at  RF). 

The  researches  mentioned  so  far  primarily  performed  interference  suppression  in  the  time- 
domain  using  the  estimation- type  predictive  filter.  However,  throughout  the  1980’s,  work  has  been 
accomplished  in  the  frequency  domain  using  transform-domain  interference  excisors. 

In  1980,  Milstein  and  Das  provided  a  detailed  analysis  of  the  performance  of  systems  in  which 
the  received  signal  is  Fourier  transformed  in  real  time  (usually  with  a  surface  acoustic  wave  device) 
and  then  filtered  by  a  multiplication  of  the  transformed  signal  by  an  appropriate  transfer  function 
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(Milstein  and  Das,  1980).  In  1989,  Gevargiz  ,  Das  and  Milstein  developed  the  transform  domain 
processing  and  presented  two  transfer  domain  processing  techniques.  In  the  first  technique,  the 
narrow  band  interference  is  detected  and  excised  in  the  transform  domain  by  using  an  adaptive 
notch  filter.  In  the  second  technique,  the  interference  is  suppressed  using  soft-limiting  in  the  trans¬ 
form  domain  (Gevargiz  et  al,  1989). 

Usage  of  a  Hamming  window  was  proposed  to  concentrate  the  energy  of  the  narrow-band 
interference  in  a  smaller  fraction  of  the  spectrum  by  Davidovici  and  Kanterakis  (Davidovici  and 
Kanterakis,  1989).  The  Hamming  Window  effectively  reduced  the  interferences  sidelobe  at  the 
expense  of  a  wider  main  lobe.  They  presented  the  overlap  and  save  algorithm  in  order  to  eliminate 
transient  effects  and  inter-symbol  interference.  Due  to  the  overlap  and  save  method  the  system’s 
filter  complexity  is  double  that  of  Milstein ’s. 

Saulnier  suggested  the  use  of  transform-domain  adaptive  filters  which  combined  some  features 
of  the  time-domain  adaptive  filter  and  the  transform  domain  excisor  (Saulnier,  1992).  Weight  leak¬ 
age  is  employed  to  allow  jammer  suppression  while  preserving  the  desired  DS  signal,  which  made 
the  desired  signal  power  lower  to  indicate  the  high-power  jammer.  Bit  error  rate  (BER)  results 
obtained  by  computer  simulation  were  presented  to  illustrate  performance  for  a  single-tone  jammer 
and  for  a  jammer  consisting  of  a  second  DS  signal  having  a  lower  chip  rate.  These  results  were 
compared  to  those  for  a  transform  domain  excisor.  Figure  2.5  show  the  simulation  block  diagram 
for  comparison  the  adaptive  filter  with  excisor  in  transform  domain. 

Another  paper  by  Reed  and  Feintuch  reported  a  comparison  between  implementations  in 
the  time  and  frequency  domains  (Reed  and  Feintuch,  1981).  In  addition,  a  thorough  review  paper 
covers  many  of  the  areas  of  interest  within  the  interference  rejection  domain  of  spread  spectrum 
communication  (Milstein,  1988) 

Relevant  research  at  the  Air  Force  Institute  of  Technology  (AFIT)  includes  the  following 
work.  Shepard’s  1982  master’s  thesis  evaluated  the  performance  of  a  DS  (Direct  Sequence)  spread 
spectrum  receiver  preceded  by  an  adaptive  interference  suppression  (AIS)  filter  (Shephard,  1982). 
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Figure  2.5  Simulation  block  diagram 


In  1984,  Way  developed  a  simulation  program  to  implement  a  19-tap  interference  suppression  filter 
for  DS-SS  signals.  Way  implemented  a  different  adaptive  algorithm:  a  soft- constraint  version  of 
the  least-mean-square  (LMS)  algorithm.  He  examined  the  filter’s  effectiveness  versus  single-tone 
jammers  and  frequency  hopping  tone  jammers.  He  likewise  evaluated  the  performance  of  the  filter 
via  SNR  improvement  (Way,  1984).  In  1990,  Mikulanicz  provided  a  performance  evaluation  of 
a  128-tap  PTF  (Programmed  Transversal  Filter).  He  demonstrated  the  effectiveness  in  creating 
assorted  FIR  bandpass  filter  responses  using  the  rectangular,  Hamming  ,  and  Hanning  window 
distribution  (Mikulanicz,  1990). 

2.3.2  Frequency  Hopping  Spread  Spectrum  with  Adaptive  Beamforming.  Adaptive  beam¬ 
forming  and  frequency  hopping  (FH)  over  a  wide  band  are  two  of  the  most  powerful  methods  for 
interference  rejection.  The  application  of  FH  modulation  in  an  antenna  array  will  improve  sub¬ 
stantially  the  SINR  (Signal  to  Interference  plus  Noise  Ratio)  as  a  consequence  of  the  increase  of 
the  resolution  and  the  interference  rejection. 

Compton  studied  the  adverse  effects  of  FH  modulation  in  an  adaptive  array  based  on  the 
LMS  algorithm  (Compton,  1985).  Bakhru  proposed  a  specific  method  for  adaptive  arrays  using 
FH  signals,  so  called  maximum  algorithm.  It  is  an  adaptive  algorithm  that  suppresses  interference 
before  it  enters  the  receivers  of  a  frequency  hopping  communication  system.  Thus,  the  algorithm 
provides  a  processing  gain  that  supplements  the  inherent  processing  gain  of  the  frequency  hopping 
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system.  The  algorithm  discriminates  between  the  desired  signal  and  the  interference  on  the  basis 
of  the  distinct  spectral  characteristics  of  frequency  hopping  signals.  The  maximum  algorithm  is  so 
named  because  the  desired  signal  is  enhanced  and  the  interference  is  suppressed  simultaneously  by 
parallel  sub-processors.  The  algorithm  is  blind  in  the  sense  that  neither  the  direction  of  the  desired 
signal  nor  its  waveform  needs  to  be  known  (Bakhru  and  Torrieri,  1984).  Nonetheless  any  adaptive 
algorithm  presents  some  discontinuities  when  used  with  FH  modulated  signals.  The  reason  is  that 
the  changes  in  the  signal  frequency  due  to  FH  are  seen  by  the  algorithm  as  changes  in  the  direction 
of  arrival. 

Torrieri  suggested  three  different  techniques  of  frequency  compensation  for  the  maximum  al¬ 
gorithm  to  solve  this  problem  (Torrieri  and  Bakhru,  1987).  The  parameter  dependent  processing 
uses  an  adaptive  filter  behind  each  antenna  element.  Each  adaptive  filter  has  enough  adjustable 
parameters  to  allow  the  formation  of  nulls  in  the  directions  of  the  sources  of  interference  for  all 
frequencies.  This  processing  is  the  most  complicated  to  implement  but  presents  faster  convergence. 
Spectral  processing  is  based  upon  dividing  that  total  hopping  band  into  a  number  of  spectral  re¬ 
gions  called  bins  and  adapting  the  weight  vector  independently  each  time  the  carrier  frequency  of 
the  frequency  hopping  signal  is  in  one  of  the  bins.  This  processing  is  the  simplest  to  implement  but 
the  achieved  improvement  is  not  significant.  An  anticipative  adaptive  system  begins  adaptation 
toward  the  optimal  weights  for  a  carrier  frequency  before  that  frequency  is  transmitted.  A  time 
advanced  frequency  hopping  replica  hops  approximately  one  hop  duration  ahead  of  the  replica  hops 
for  dehopping  the  received  signal.  While  the  main  adaptive  filter  produces  the  output,  the  auxiliary 
filter  adapts  its  weights  to  suppress  the  next  carrier  frequency.  After  each  hop,  the  weight  values 
associated  with  the  new  carrier  frequency  are  transferred  from  the  auxiliary  filter  to  the  main  filter. 
Anticipate  processing  provides  the  fastest  convergence  to  the  steady  state,  but  exhibits  the  largest 
variation  in  the  steady  state  SINK. 

Najar  and  Lagunas  further  developed  the  anticipative  processing  by  adding  a  generalized 
sidelobe  canceller  (GSC).  His  approach  to  the  optimum  solution  consists  two  different  stages.  The 
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first  stage,  named  the  anticipative  stage,  is  devoted  to  cancel  interference  at  fixed  frequencies  that 
are  already  present  or  active  at  the  frequency  of  interest  at  the  hop  time.  The  second  state,  named 
GSC  stage,  is  used  for  combating  interferences  that  are  not  present  at  the  frequency  of  interest  at 
the  time  of  the  frequency  hop  (Najar  and  Lagunas,  1995).  They  showed  the  grating  lobe  is  reduced 
significantly  in  the  mean  array  factor  when  the  FH  signals  is  used  with  antenna  array. 

2.4  Conclusion 

This  chapter  presented  a  brief  historical  review  of  developments  in  three  areas  that  are  closely 
related  in  so  far  as  the  subject  matter  of  this  research.  Those  areas  are  adaptive  filtering  algorithm, 
adaptive  beamforming  for  wideband  interference  suppression.  Many  researchers  have  investigated 
the  adaptive  filters  for  narrowband  interference  suppression  and  adaptive  beamforming  for  wide¬ 
band  interference  suppression.  Special  attention  was  devoted  to  adaptive  filters  in  the  DS-SS 
system.  In  the  FH-SS  signal  environments,  the  researches  relating  adaptive  antenna  array  and 
FH  modulation  for  interference  suppression  have  been  developed  because  the  antenna  array  and 
FH  modulation  is  closely  related  with  each  other.  In  the  next  chapter,  we  derive  the  theories  of 
adaptive  filter  and  beamforming. 
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III.  Adaptive  Signal  Processing 


3.1  Background 

3.1.1  Overview.  This  chapter  discusses  the  similarity  of  temporal  and  spatial  filters. 
Then,  the  time  domain  adaptive  filters  are  considered  using  LMS  adaptive  algorithm  and  the  spatial 
filter  (i.e.,  adaptive  beamforming)  is  discussed  using  the  linear  constrained  minimum  variance 
(LCMV)  algorithm  with  generalized  sidelobe  canceller  (GSC).  An  optimum  value  and  LMS  estimate 
of  weights  are  also  computed.  Finally,  the  adaptive  beamformer  in  the  presence  of  correlated  signals 
is  considered. 

3.1.2  Notation.  Since  vectors  and  matrices  are  used  throughout  this  research,  it  is 
important  to  establish  the  notation  that  will  be  used.  Vectors  are  always  represented  with  lowercase 
boldface  symbols,  for  example,  w,  and  are  assumed  to  be  column  vector.  Matrices  are  denoted 
by  boldface  upper  case  symbols,  for  example  C.  Superscripts  *,  T,  if,  and  -1  represent  complex 
conjugate,  transpose,  complex  conjugate  transpose,  and  matrix  inverse,  respectively. 

3.1.3  Definition,  Signals  may  be  classified  as  either  narrowband  or  wideband.  Narrow- 
band  is  defined  in  terms  of  the  fractional  bandwidth  of.  the  signal.  The  fractional  bandwidth  is 
the  signal  bandwidth  as  a  percentage  of  the  carrier  frequency.  Signals  whose  fractional  bandwidths 
are  much  less  than  2  percent  will  be  characterized  as  narrowband,  while  those  with  fractional 
bandwidths  much  greater  than  2  percent  will  be  called  wideband  (Haykin  and  Steinhardt,  1992). 

3.1.4  Temporal  Filter  versus  Spatial  Filter.  Although  the  temporal  filter  problem  in 
additive  receiver  noise  and  the  spatial  filter  problem  corrupted  by  additive  sensor  noise  arise  in 
different  application  areas,  their  mathematical  formulations  and  procedures  for  their  solution  are 
indeed  similar. 

For  the  temporal  filtering  of  data,  we  propose  using  a  transversal  filter  of  length  M  as  indicated 
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in  Fig.  2.1.  The  filter  output,  in  response  to  the  tap  inputs,  is  given  by 

M-l 

fc=0 

where  Wk  is  A:th  tap  weight  and  14(71)  is  tap  input  in  time  n.  For  the  special  case  of  a  sinusoidal 
excitation 

u{n)  =  (3.2) 


we  may  rewrite  Eq.  3.1  as 


Af-l 


y{n)  =  ^ 

fc=0 


(3.3) 


where  u  is  the  angular  frequency  of  the  excitation,  which  is  normalized  with  respect  to  the  sampling 
rate. 


Consider  next  the  spatial  analog  of  this  temporal  problem.  Fig.  3.1  depicts  a  receiving  linear 
array  of  uniformly  spaced  sensors  labeled  1,-  ■  ■,  M.  Let  d  denote  the  separation  between  adjacent 
elements  of  the  array.  In  Fig.  3.1  a  single  plane  wave  was  impinging  on  the  array  at  angle  of 
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incidence  9,  measured  with  respect  to  the  boresight.  It  is  assumed  that  the  inter-element  distance 
is  a  half  wavelength  so  as  to  avoid  the  appearance  of  grating  lobes  (Skolnik,  1980).  The  resulting 
beamforming  output  is  given  by 


M-l 

y{n)  =  u{n)  ^  (3.4) 

k=0 

where  the  direction  of  arrival  is  defined  by  the  electrical  angle  4)  =  sin  6  that  is  related  to  the 
angle  of  incidence  9  and  the  antenna  interelements  distance  d,  u{n)  is  the  electrical  signal  picked  up 
by  the  antenna  element  labeled  ui  in  Fig.  3.1  that  is  treated  as  the  point  of  reference,  and  the  Wk 
denote  the  element  weights  of  the  beamformer.  The  important  point  to  note  is  the  mathematical 
similarity  between  the  temporal  model  of  Eq.  3.3  and  the  spatial  model  of  Eq.  3.4. 

3.2  Adaptive  Filter  :  Temporal  Filter 

An  adaptive  filter  converges  to  the  optimum  Wiener  solution  in  a  stationary  environment. 
The  LMS  algorithm  starts  from  some  predetermined  set  of  initial  conditions  and  is  self-designing. 
It  relies  for  its  operation  on  a  recursive  algorithm,  which  makes  it  possible  for  the  filter  to  perform 
satisfactorily  in  an  environment  where  complete  knowledge  of  the  relevant  signal  characteristics  is 
not  available.  On  the  other  hand,  the  design  of  a  Wiener  filter  requires  prior  information  about 
the  statistics  of  the  data  to  be  processed  (Widrow,  1985).  For  real-time  operation,  this  procedure 
has  the  disadvantage  of  requiring  excessively  elaborate  and  costly  hardware. 

3.2.1  Input  Signal  and  Weight  Vector.  The  operation  of  a  linear  adaptive  filtering  algo¬ 
rithm  involves  two  basic  processes:  (1)  a,  filtering  process  designed  to  produce  an  output  in  response 
to  sequence  of  input  data,  and  (2)  an  adaptive  process.  The  transversal  filter  type  is  used,  which 
referred  to  as  a  tapped-delay  line  filter,  as  depicted  in  Fig.  2.1.  For  the  transversal  filter,  we  obtain 
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the  input-output  relationships  as  follows: 

Af-l 

2/ W  =  52  “  ^)-  (3-^) 

k=0 

In  vector  notation,  this  can  be  written  as 

y{n)  =  w"u(n),  (3.6) 

where  the  input  and  weight  vector  are  expressed  by  vector  notation  as, 

u(n)  =  [tt(n),M(n —  !),•••  ,M(n  —  M -|- 1)]^ 
w  =  [wi,W2,---  ,wm]'^- 

3.2.2  Minimum  Mean  Squared  Error.  For  the  case  of  stationary  inputs,  the  performance 
function  is  defined  as  the  mean  square  error  which  is  the  mean  squared  value  of  the  difference 
between  the  desired  response  and  the  transversal  filter  output.  As  in  Fig.  2.2  the  error  signal  is 

e(n)  =  d{n)  -  y{n) 

=  d{n)  —  wr^u. 

Squaring  the  error  signal  to  obtain  the  squared  error,  one  obtains 

e(n)^  =  d{nY  -f  w^uu^w  —  2d{n)u^w.  (3.7) 

Then  taking  the  expected  value,  the  MSB  (Mean  Squared  Error)  becomes 

MSE  =  B{d(ra)2}-fw^B{uu^}w-2B{d(ri,)u^}w 
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£{d^}  +  w^Rw  -  2pw, 


(3.8) 


where  R  is  defined  as  the  square  matrix  : 

R  =  i7{uu^}.  (3.9) 

Correspondingly,  p  is  denoted  by  the  M-by-1  cross  correlation  vector  between  the  tap  inputs  of 
the  filter  and  the  desired  response  d{n)  : 

p  =  £l{ud(n)*}.  (3.10) 

Many  useful  adaptive  processes  work  by  seeking  the  weight  vector  which  minimizes  the  perfor¬ 
mance  function  using  gradient  search  methods.  The  gradient  of  the  mean  squared  error  performance 
surface,  designated  V,  can  be  obtained  by  dififerentiating  Eq.  3.8  to  obtain  the  column  vector  : 

V=  2Rw-2p. 

3.2.3  LMS  Adaptation  Algorithm.  Since  the  gradient  vector  must  be  estimated  from  the 
available  data,  the  simplest  estimate  of  gradient  vector  is  the  instantaneous  estimates  based  on 
sample  values  of  the  tap-input  vector  and  desired  response  : 

V  =  2Rw  —  2p  =  2u(ri)u'^(n)w  —  2vi{n)d*[n) 

Let  w(n)  denote  the  value  of  the  tap  weight  vector  at  time  n.  The  updated  value  of  the 
tap-weight  vector  at  time  n  -f  1  is  computed  by  using  the  simple  recursive  relation 

■w{n  -k  1)  =  w{n)  -k  fiii{n)\d* {n)  —  u^(n)w]  (3-11) 
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3.2-4  Summary.  Adaptive  temporal  filter  have  a  built  in  mechanism  for  the  automatic 
adjustment  of  its  weight  in  response  to  statistical  variations  of  the  environment  in  which  the  filter 
operates.  The  weight  adjustments  are  made  iteratively  by  using  standard  adaptive  algorithms  that 
seeks  an  optimum  value.  The  filter  input  consisting  of  a  vector  of  uniformly  spaced  samples  taken 
from  a  long  data  stream  is  used  to  determine  weights  that  minimize  the  performance  function 
calculated  by  gradient  methods.  The  updated  weight  vector  is  computed  by  using  the  simple 
recursive  relation  as  depicted  in  Eq.  3.11. 

3.3  Adaptive  Beamforming 

3.3.1  Introduction.  A  beamformer  is  a  signal  processor  used  in  conjunction  with  a  set  of 
antennas  that  are  spatially  separated.  The  beamformer  output  is  simply  a  weighted  combination 
of  the  outputs  of  the  set  of  antennas.  Usually  the  goal  of  beamforming  is  spatial  filtering,  that  is, 
separation  of  signals  which  have  similar  temporal  frequency  content  but  originate  from  different 
spatial  locations.  Here,  we  use  digital  beamforming  that  detects  and  digitizes  the  received  signal 
at  the  element  level  via  discrete  processing  techniques. 

In  an  adaptive  beamformer  the  weights  are  adjusted  or  adapted  in  response  to  the  data  re¬ 
ceived  at  the  antennas  to  optimize  the  beamformer’s  spatial  response.  Changing  the  beamformer’s 
spatial  response  is  as  simple  as  changing  the  weights.  The  weights  in  an  adaptive  beamformer  are 
controlled  by  an  adaptive  algorithm.  All  adaptive  antennas  to  date  have  been  arrays,  because  the 
pattern  of  an  array  is  easily  controlled  by  adjusting  the  amplitude  and  phase  of  the  signal  from  each 
element  before  combining  the  signals.  Adaptive  antennas  are  useful  in  radar  and  communication 
systems  that  are  subject  to  interference  and  jamming.  They  change  their  patterns  in  a  way  that 
optimizes  the  signal-to-interference-plus-noise  ratio  at  the  array  output. 
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3.3.2  Input  Signal  and  Weight  Vector  .  Fig.  1.1  illustrated  a  beamformer  typically 
used  for  processing  narrowband  signals.  The  antenna  outputs  are  sampled  and  are  represented 
in  discrete  time  for  notational  convenience  and  consistency.  The  M  antenna  sensor  outputs  are 
weighted  and  then  summed  to  compute  the  beamformer  output.  This  method  is  called  as  delay- 
and-sum  beamformer.  Each  sensor  is  assumed  to  have  any  necessary  receiver  electronics  and  an 
A/D  converter.  The  discrete  time  delay  and  sum  beamformer  2/(n)  is  given  by 

M 

yi^)  =  H  'W^Umin),  (3.12) 

m=l 

where  Wm  and  Um{n)  are  the  weight  and  input  data,  respectively,  in  the  mth  antenna  channel. 

A  beamforming  structure  commonly  utilized  to  process  wideband  signals  is  depicted  in  Fig.  2.4. 
In  this  case  there  are  tap  delay  lines  in  each  channel.  The  output,  y{n)  is  expressed  as 

M-lL-l 

2/W  =  S  (3.13) 

m=l  i=0 

where  L  is  the  number  of  taps  in  each  of  the  M  channels  and  Wm,i  is  the  weight  applied  to  the 
Zth  tap  of  the  mth  channel.  This  means  that  the  output  of  each  sensor  is  then  passed  through  a 
FIR  filter  having  L  weights.  The  time  domain  filtering  is  intended  to  provide  some  rejection  of 
interference  not  lying  in  the  proper  temporal  frequency  region.  The  resulting  FIR  filter  describing 
the  relation  between  array  input  and  output  has  tap  values  equal  to  the  sum  of  component-filter  tap 
values  occurring  at  the  temporal  delay.  Once  each  antenna  sensor’s  output  delayed  by  an  amount 
appropriate  for  the  propagation  direction  is  passed  through  a  FIR  filter  whose  weights  are  selected 
individually. 

Both  Eq.  3.12  and  Eq.  3.13  are  compactly  written  as  the  inner  product 

y{n)  =  w^u{n),  (3.14) 
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where  w  is  the  vector  of  beamformer  weights  and  u(n)  is  vector  of  data. 

In  vector  notation,  the  complex  conjugate  transpose  of  the  iV-dimensional  vector  of  signals 
contained  in  the  beamformer  at  the  nth  time  sample  is  given  by 

u^(n)  =  [uf  (n),  uf(n),  uf_i(n)].  (3.15) 

Let  the  dimension  of  w  and  u(n)  be  N  where  N  =  M  in  the  narrowband  case  (Eq.  3.12)  and 
N  =  ML  in  the  wideband  case  (Eq.  3.13).  Note  that  the  form  of  both  equation  implies  that  the 
complex  conjugate  of  the  weights  are  actually  applied  to  the  data  samples  within  the  beamformer. 


3.3.3  Antenna  Array  Response  Vector.  In  the  narrowband  case,  the  frequency  response 
of  an  antenna  array  with  antenna  element  weights  Wm  and  a  propagation  delay  Tp  seconds  is  given 
by, 

r  (wc,  Op)  =  w*d(wc,  Op),  (3.16) 

where 

w  =  [wi,  W2,  •••,  wm]^. 

d{u„0p)  =  [1 

r(u)c,0p)  represents  the  response  of  the  antenna  to  frequency  and  d{u)c,0p)  is  a  vector  describing 
the  phase  at  each  element  in  the  antenna  array  relative  to  the  element  associated  with  Wi.  Tp  is 
^sin(0p),  where  c  is  propagation  velocity  and  d  is  antenna  interelement  space  and  Op  is  incidental 
direction  for  signals. 

Assume  that  the  array  spatially  samples  a  propagating  wave  of  frequency  Wc  in  the  direction 
Op.  For  convenience  let  the  phase  due  to  propagation  be  referenced  to  zero  at  the  first  sensor  and 
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let  ai(0p,uJc)  be  the  response  of  the  first  antenna  element  as  a  function  of  direction  and  frequency. 
Assuming  identical  antennas, 


M-l 

y{n)  = 

m=0 

=  u{n)r(6p,(Jc)- 

The  term  r[6p,(jJc)  is  the  beamformer  response  and  the  elements  of  d(&p,a;c)  correspond  to  the 
complex  exponentials  1  <  m  <  M.  The  vector  d(^p,Wc)  is  termed  the 

array  response  vector.  Assuming  ai(^p,Wc)=l,  d(6p,(Jc)  is  written  as 

d(^p,a;c)  =  [l  •••  gAM-iVcTpjit,  (3,17) 

In  this  case  we  identify  the  relationship  between  temporal  frequency  u  in  d(u;)  (FIR  filter)  and 
direction  6p  (beamformer)  as  w  =  Wcf  sin  (Op). 

3.3.4  Basic  Concept  of  LCMV  Beamforming  . 

3. 3. 4-1  Linear  Constraints.  We  will  use  the  constraint,  in  which  the  beamformer 
response  at  a  specified  frequency  and  direction  of  desired  signal  must  meet  some  gain  and  phase 
requirements.  The  response  is  constrained  at  radian  frequency  0)^  and  angle  Op  as 

w"d(0p,a;c)  = /i,  (3.18) 

where  /i  is  the  desired  complex  response.  Several  linear  constraints  can  be  expressed  as, 

C^w  =  f,  (3.19) 
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where  C  is  the  constraint  matrix  and  f  is  the  response  vector. 

In  the  narrowband  case,  the  matrix  C  is  chosen  to  contain  r  linear  constraints  so  the  above 
equation  represents  r  linearly  independent  columns  in  M  unknowns.  The  weight  vector  w  has  length 
M,  while  C  is  an  M  xr  matrix.  If  M  =  r,  then  w  is  uniquely  determined  by  the  constraints.  To 
ensure  there  is  a  w  which  satisfies  the  constraints,  r  is  chosen  to  be  less  than  M.  Here  we  consider 
the  desired  signal  gain-only  constraint, 

C  =  d(0p,n;e)  (3.20) 

where 

d(dp,u;e)  =  •••  ^  ^3_2i) 

d(^p,a;c)  is  M  X  1  and  Tp  is  ^  sin^p,  c  is  propagation  velocity  and  d  is  antenna  interelement  space 
and  9p  is  incidental  direction  for  signals. 

In  the  wideband  case,  for  steer-direction  gain-only  constraints, 


f 

=  [l,0...,0f 

(3.22) 

d(flp,  cUc) 

Om 

Om 

Om 

d(0p,  u)c)  ... 

Om 

(3.23) 

Om 

Om 

d(^p,  ^ 

The  d(0p,Wc)  is  ilf  X  1,  and  the  column  vector  Om  contains  M  zeros.  The  T-dimensional  gain 
vector  f  describes  the  frequency  response  of  the  beamformer  to  a  signal  impinging  on  the  array 
from  the  desired  look  direction  (Frost,  1972).  The  constraint  matrix  C  is  ML  x  L  dimensional.  If 
f  is  a  vector  which  contains  a  single  one  and  L  —  1  zeros,  then  the  temporal  response  will  be  flat 
and  the  system  will  not  distort  any  signal  incident  on  the  array  from  the  look  direction.  Op. 
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3. 3. 4-2  LCMV  Approach.  Variance  minimization  effectively  minimizes  the  interfer¬ 
ence  and  noise  power  at  the  beamformer  output  while  the  constraints  preserve  the  signal  of  interest. 
Formally  stated,  we  derive 

min(Po)  =  min(w^R„w),  subject  to  C^w  =  f,  (3-24) 

where  Pg  represents  the  output  power.  Ru  =  P{uu®}  represents  the  exact  covariance  matrix  of 
the  received  signal  u.  The  Eq.  3.24  is  to  find  out  the  weight  vector  for  minimum  output  Pg.  The 
solution  is  derived  in  Appendix  A  using  Lagrange  multipliers,  given  by, 

wopt  =  R-'C[C®R-'C]-if.  (3.25) 

The  inverse  exists  because  C  is  full  rank  and  R^  is  positive  definite  because  the  data  always  contains 
an  uncorrelated  noise  component. 

3.3.5  L  CMV-  GSC  Beamformer  . 

3. 3. 5.1  Basic  Concept  of  GSC.  Fig.  3.2  shows  the  generalized  sidelobe  canceler 
(GSC)  block  diagram.  Note  that  the  input  vector  is  the  M-dimensional  stacked  data  vector  u(n) 
for  the  narrowband  signal.  For  steer  direction  gain  only  constraints  (r  =  1),  the  linear  constraint 
is  defined  in  Eq.  3.20.  Let  the  columns  of  an  M  x  (M  —  r)  matrix  Cg  be  defined  as  a  basis  for 
the  orthogonal  complement  of  the  space  spanned  by  the  columns  of  matrix  C  which  is  composed 
of  orthonormal  columns.  Using  the  definition  of  an  orthogonal  complement,  we  may  thus  write 

C^Cg  =  0.  (3.26) 
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e(n) 


Figure  3.2  Block  Diagram  for  a  Generalized  Side-Lobe  Canceler  (GSC) 

Define  the  M-hy-M  partitioned  matrix 

U  =  [C  :  Ca],  (3.27) 

whose  columns  span  the  entire  M-dimensional  signal  space. 

Let  the  v  be  partitioned  in  compatible  to  Eq.  3.27, 

q  ■ 

...  .  (3.28) 

-w„. 

where  q  is  r-by-1  vector,  and  the  (M  —  7-)-by-l  vector  Wo  is  the  portion  of  the  weight  vector  w 
that  is  not  affected  by  the  constraints. 

So  that  the  M-by-1  weight  vector  is  represented  by  U  and  v  where 

w  =  Uv,  (3.29) 
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and  U  ^  exists  because  the  determinant  U  is  not  zero, 


V  =  U-^w.  (3.30) 

We  may  use  Eq.  3.27  and  Eq.  3.28  in  Eq.  3.29  to  obtain 

■  q  ■ 

w  =  [C:C„]  ...  (3.31) 

.-Wa. 

=  Cq-C„w„.  (3.32) 

Now,  inserting  Eq.  3.32  into  Eq.  3.19,  we  obtain 

C^Cq  -  C^CaWa  =  f.  (3.33) 

Through  Eq.  3.26,  the  Eq.  3.33  will  be, 

C^Cq  =  f.  (3.34) 

Therefore,  solving  for  the  vector  q, 

q=(C^C)-if,  (3.35) 

which  shows  the  linear  constraint  doesn’t  affect  Wa. 

Define  the  non-adaptive  beamformer  component  represented  by,  which  is  the  first  part  of 
Eq.  3.32, 

w,  =  Cq  =  (3.36) 
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Substituting  Eq.  3.36  into  Eq.  3.32  yields  the  follow  equation  relating, 


w  =  Wg  -  C„Wa. 


(3.37) 


Substituting  Eq.  3.37  into  Eq.  3.19, 


C^Wg  -  C^CaVfa  =  f.  (3.38) 

By  virtue  of  Eq.  3.26, 

C^Wg  =  f,  (3.39) 

which  shows  the  non-adaptive  component  Wg  of  the  weights  w  satisfies  the  linear  constraints.  The 
weight  vector  Wg  is  termed  the  quiescent  weight  vector  because  it  is  optimum  LCMV  weight  vector 
when  the  environment  is  quiet.  As  the  results,  the  GSC  has  only  one  data  dependent  element  Wg 
due  to  the  fact  that  both  the  quiescent  beamformer  Wg  in  the  upper  path  and  the  blocking  matrix 
Cg  in  the  lower  path  depend  only  on  the  constraint  equations. 

Note  that  the  constraint  does  not  affect  Wg  because  Wg  embodies  the  available  degrees  of 
freedom  in  w.  The  second  term  on  the  right  side  of  Eq.  3.37  is  incorporated  in  the  lower  branch 
of  the  GSC  and  is  essentially  an  unconstrained  adaptive  beamforming  ,  which  have  the  M  —  r  the 
number  of  adaptive  weights.  The  constraint  matrix  C  preserves  the  desired  signal.  But  the  matrix 
Cg  blocks  any  portion  of  the  data  that  is  constrained  in  the  space  spanned  by  the  columns  of  C. 
Consequently,  no  part  of  the  received  signal  subspace  is  allowed  to  pass  through  the  lower  branch 
of  the  GSC.  The  matrix  Cg  is  for  this  reason  called  the  blocking  matrix. 

3. 3. 5. 2  LMS  Estimate  of  Weights  for  LCMV-GSC  Beamforming.  We  perform  the 
unconstrained  minimization  of  the  mean-squared  value  to  minimize  the  effect  of  the  jammer  by 
adjusting  the  weight  vector  Wg  using  LMS  algorithm.  According  to  Eq.  3.14,  the  beamformer 
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output  y{n)  is  expressed  by  the  inner  product, 

y{n)  =  ■vv®u(n),  (3.40) 

where  u[n)  is  the  electrical  signal  picked  up  by  first  antenna  element  of  the  linear  array  in  Fig.  1.1 
and  Fig.  2.4  at  time  n.  By  substituting  Eq.  3.37  in  Eq.  3.40,  the  beamformer  output  in  Fig.  3.2  is, 

y{n)  =  wf  u(n)  -  wf  Cf  u(n.).  (3.41) 

In  Eq.  3.41,  the  inner  product  w^u(n)  plays  the  role  of  quiescent  response: 

q{n)  =  wf  u(n).  (3.42) 

Similary,  the  matrix  product  C^u(ra)  plays  the  role  of  input  vector  for  the  adjustable  weight  vector 
Wfl.  To  emphasize  this  point,  we  let 

x(n)  =  C®u(n).  (3.43) 

We  are  now  ready  to  formulate  the  LMS  algorithm  for  the  adaption  of  weight  vector  Wa(n) 
in  the  GSC.  The  derivation  is  based  on  Section  3.2,  especially  Eq.  3.11.  We  may  write 

Wa(n+1)  =  Wa(n) +  /tCfu(n)e*(n) 

=  Wa(n)  +  /tCf  u(n)(wf  u(n)  -  wf  Cf  u(n))* 

=  Wa(ri) +  /xCfu(w)u(n)^(wg  -  CaWo(n)).  (3-44) 

3.3.5.S  Optimum  Value  of  Weights  for  LCMV-GSC  Beamforming.  Based  upon 
this  explanation  of  the  GSC,  we  can  now  express  the  LCMV  beamformer  with  the  following 
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unconstrained  minimization  equation  (Haykin  and  Steinhardt,  1992), 


Pa  -  min(W5  -  CaWa)-^R„(Wg  -  CaWfl), 


(3.45) 


or 


Pa  =  min(wf  RuWg  -  w^R„CaW«  -  wf  Cf  R^w^  +  wf  Cf  R„CaWa).  (3.46) 


Completing  the  square  we  have 


Pa  =  min  [w«  -  (Cf  R„C«)-iCf  R„w,]Cf  R„C„[w,  -  (Cf R„C„)-iCfR„w,] 

+wf  R„w,  -  wf  R„C„(Cf  R„Ca)-iCf  R„Wj.  (3.47) 

The  matrix  C„  is  full  rank  and  R„  is  positive  definite  so  C^RuC„  is  positive  definite  and  the 
inverse  exists.  Eq.  3.47  is  minimized  when 

w„  =  (Cf  R„C„)-iCf  R„w,.  (3.48) 

The  minimum  output  power  is  obtained  using  this  set  of  weights  and  is  given  by 

Po  =  w"R„w,  -  wf  R„Ca(Cf  R„Ca)“^C^R„Wj.  (3.49) 


Suppose  u  =  s  +  j  where  s  is  the  component  of  the  data  due  to  the  desired  signal  and  j  is  the 
component  due  to  jammer  signal.  Assume  that  the  constraints  are  chosen  to  preserve  the  desired 
signal,  that  is,  s  lies  in  the  space  spanned  by  the  columns  of  C  so  Cf^Rs  =  0.  This  implies  that 
the  beamformer  output  is 

y  =  w^s  +  (wg  -  C„Wa)^j.  (3.50) 
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If  we  further  assume  that  the  signal  and  interference  are  uncorrelated,  then  the  output  power  or 
variance  is  expressed  as 

E{\y\'^}  =  -wfRsWq  +  w^RjW,  (3.51) 

where  w^RjW  is  the  jamming  signal  component  of  the  output  power.  Therefore,  minimizing  the 
total  output  power  subject  to  C'^w  =  f  is  equivalent  to  minimizing  the  jammer  output  power 
subject  to  the  same  constraints.  Substitute  the  GSC  representation  for  w  into  Eq.  3.51  with  Wa 
given  by  Eq.  3.48  and  utilize  the  identity  C^Rg  =  0  to  express  the  minimum  output  power  Po  as 
a  sum  of  signal  output  power,  Pg,  and  jammer  output  power,  Pj  where 

Pg  =  wfRgW, 

P„  =  w^RjW.  (3.52) 

The  mean  squared  error  (MSE)  between  the  desired  signal  and  the  beamformer  output  is  easily 
derived  using  the  GSC  representation.  We  assume  the  constraints  are  chosen  so  that  the  beamformer 
output  in  the  absence  of  interference  and  noise  is  equal  to  the  desired  signal,  yd-  that  is,  yd  =  w^s. 
Defining  the  MSE  as 

MSE  =  E{\yd-y\^),  (3.53) 

and  substituting  Eq.  3.50  into  Eq.  3.53  yields 

MSE  =  P{|w"jp} 

=  w^R^w.  (3.54) 

Comparison  of  Eq.  3.52  and  Eq.  3.54  indicates  that  the  LCMV  criterion  is  equivalent  to  a  minimum 
MSE  criterion. 
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Examination  of  Eq.  3.48  reveals  the  similarity  with  the  standard  Wiener  solution.  Consider 


the  auto-covariance  of  x(n)  in  Fig.  3.2, 


Rx  =  E{xx^}  =  Cf  R^Ca, 


(3.55) 


because  x(n)  =  Cfu(  n). 

Now,  consider  the  cross  covariance  of  the  signals  x(n)  and  q{n), 

(3.56) 

(3.57) 

(3.58) 


rx,  =  E{x{n)q{n)^} 

=  •®{Cfu(n)u^(n)w,} 

=  CfR„w,. 


Comparing  Eq.  3.55  and  Eq.  3.58  with  Eq.  3.48  ,  we  can  rewrite  Eq.  3.48  as, 

Wa  =  Rx ’■R«rxg,  (3.59) 

which  shows  the  similarity  between  the  standard  Wiener  solution  and  the  solution  implemented 
with  a  GSC. 

3.3.6  Summary.  A  beamformer  forms  a  weighted  combination  of  the  outputs  of  spatially 
separate  antennas  in  order  to  separate  propagating  signals  that  originate  from  different  locations. 
The  weights  determine  the  spatial  and  temporal  filtering  characteristics  of  the  beamformer.  In 
an  adaptive  beamformer  the  weights  are  adjusted  in  response  to  the  data  received  at  the  antenna 
for  the  purpose  of  optimizing  the  beamformer’s  response.  A  common  criterion  for  optimizing  the 
weights  is  minimization  of  beamformer  output  power  or  variance.  The  goal  of  minimizing  output 
power  is  to  minimize  the  contributions  of  jammer  to  the  output.  The  LCMV  approach  minimizes 
output  power  subject  to  a  set  of  linear  constraints  on  the  weight  vector.  These  constraints  are 
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used  to  control  the  response  of  the  beamformer  in  the  direction  of  the  desired  signal  to  prevent  the 
weights  from  cancelling  the  desired  signal. 

The  GSC  separates  the  weight  vector  into  constrained  and  unconstrained  components.  The 
unconstrained  components  represent  the  beamformer’s  adaptive  degrees  of  freedom  and  are  adjusted 
using  standard  adaptive  algorithm.  In  order  to  suppress  the  jammer  signal  from  different  incoming 
direction,  the  optimum  weight  vector  and  LMS  estimate  of  weight  vector  were  calculated  (Eq.  3.48 
and  Eq.  3.44). 

3.4  Adaptive  Beamforming  in  the  Presence  of  Correlated  Signals 

3.4-1  Introduction.  The  behavior  of  adaptive  arrays  when  the  interference(s)  is  correlated 
with  the  desired  signal  is  of  concern  in  some  environments.  We  will  call  two  signals  fully  correlated 
if  one  is  a  scaled,  delayed  replica  of  the  other.  Correlation  can  destroy  the  performance  of  a 
constrained  adaptive  array  through  two  effects:  1)  the  beamformer  fails  to  form  deep  nulls  in 
the  directions  of  the  correlated  interferences  (Reddy  et  al,  1987)  and  2)  the  desired  signal  can 
be  partially  or  completely  cancelled  (Widrow  et  al,  1982).  This  problem  can  be  overcome  by  a 
technique  called  spatial  smoothing  which  modulates  the  interference  in  a  way  that  will  reduce  the 
correlation  with  the  desired  signal.  A  method  of  combining  spatial  smoothing  with  LMS  algorithm 
will  be  illustrated. 

3.4.2  Analysis  of  the  Decorrelation  Effect  of  Spatial  Smoothing.  In  the  GSC  depicted  in 
Fig.  3.2,  the  weights  Wq  passes  the  desired  signal  through  the  upper  branch  while  the  matrix  Ca 
blocks  it  from  the  lower  branch.  If  however  the  desired  signal  is  correlated  with  an  interference, 
the  signals  in  the  lower  branch  of  the  GSC  structure  can  be  used  to  cancel  the  desired  signal  in  the 
upper  branch.  This  is  because  the  interference  is  a  version  of  the  desired  signal  that  is  not  removed 
from  the  lower  branch  by  the  spatial  filter  Co. 
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Let  u(n)  be  the  simultaneously  sampled  vector  (snapshot)  of  array  signals. 

u(n)  =  [uf(n),  uf(n),  uf_i(n)]^.  (3.60) 

where 

(n)  =  [uixM{n),  UixM+iin),  •••,  ttixM+M-i(w)],  0  <  i  <  i  -  1.  (3.61) 

Let  the  dimension  of  u(n)  be  N  where  N  =  M  in  the  narrowband  case  (Eq.  3.12)  and  N  =  M  x  L 
in  the  wideband  case  (Eq.  3.13).  In  Eq.  3.60,  u(n)  is  iV  x  1  with  element, 

Ui(n)  =  u[nT,~7nod(i-l,M)Tp-INT{^-^)To],  l<i<N,  (3.62) 

where  mod  stands  for  modulo  and  INT  is  integer.  Ts  is  sampling  time  and  Tp  represents  the  time 
delay  due  to  propagation  and  the  delay  time  To  of  each  tap  in  the  TDL  antenna  array  is  determined 
by 

To  =  ^  (3.63) 

where  Rg  is  the  sampling  frequency  of  the  arriving  signal. 

To  simplify  our  work  we  restrict  our  analysis  to  the  case  of  two  signals.  One  of  these  is  the 
desired  signal,  and  the  other  is  a  correlated  interference.  In  Eq.  3.62,  u(n)  represents  the  desired 
signal  when  Tp  is  Td  and  interference  signal  when  Tp  is  Tj . 

The  basic  concept  of  spatial  smoothing  is  as  follows.  Assume  that  the  uniform  linear  array  of 
M  sensors  is  extended  with  additional  sensors  and  the  extended  array  is  grouped  into  sub-arrays 
of  size  M.  The  first  sub-array  is  formed  from  the  sensors  1,  •  ■  • ,  M,  and  the  second  sub-array  is 
formed  from  the  sensors  2,  •  •  • ,  M  -f- 1,  and  so  on.  Let  us  denote  the  vector  of  received  signals  at 
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the  fcth  sub- array  by  u*'.  The  feth  sub  array  and  the  i  th  element  of  the  input  vector  is, 


u\in)  =  u\nT^  —  mod{i  -b  /d  —  1,M  +  k)Tp  +INT{ 


1 


M  +  k 


)ro], 


1  <  ?■  <  iV,  (3.64) 


where  the  superscript  indicates  the  A:th  subarray  and  the  subscript  indicates  the  ith  element  of  the 
input  vector.  The  array  covariance  matrix  of  the  fcth  sub-array  is  then  given  by 


=  £{(u*=)(u*’)^} 


(3.65) 


We  define  the  spatially  smoothing  covariance  matrix  with  K  number  of  smoothing  steps  as 


1  ^ 


fc=i 


(3.66) 


Combining  Eq.  3.65  and  Eq.  3.66,  we  obtain 


R=lf;E{(u*')(u^)"} 


(3.67) 


fc=i 


where  u*’(n)  is  composed  of  desired  signal  and  interference  signal  in  the  A;th  subarray.  Hence,  we 
have 


^  +  (3.68) 

k=l 

=  i  +  ^[(u))(u^)^]  (3.69) 

fc=i 

=  Rg -b  Rj -b  R^j -f- Rjs,  (3.70) 


where  Rg  is  the  autocorrelation  of  desired  signal  and  Rj  is  the  autocorrelation  of  jamming  signal, 
and  Rgj  is  crosscorrelation  between  the  desired  signal  and  jamming  signal  in  the  A;th  sub-array. 
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We  now  study  how  progressive  spatial  smoothing  reduces  the  crosscorrelation  between  all  the 
incident  wavefronts,  and  hence  between  the  desired  signal  and  the  other  interference  wavefronts, 
and  therefore  impacts  interferences  cancellation.  If  we  assume  the  incident  wavefronts  are  frequency 
hopped  signal,  we  can  write, 

u(n)  =  u{n)e(h),  Th-i  <n<Th,  (3-71) 

where  the  signal  arrival  vector,  e(/i),  is  W  x  1  with  element  i, 

ei(h)  =  exp[-j{mod(i  -  1,  M)<f,p{h)  +  /WT(^^"(h)}],  l<i<N.  (3.72) 

Here  <j>p{h)  is  changed  to  (f>dih)  for  the  desired  signal  interelement  phase  shift  during  hop  h  and  (j>j 
for  the  interference  signal  interelement  phase  shift  during  hop  j,  i.e, 

(pdih)  =  ijJhTd 
4>j{h)  = 

where  Uh  and  ujj  is  the  desired  signal  and  jamming  signal  hopping  frequencies  respectively.  There¬ 
fore,  ed{h)  represents  the  desired  signal  and  ej(h)  does  the  jamming  signal.  Also  ^p{h)  is  the  phase 
shift  by  delay  time  To,  so. 


=  u;hTo 
=  ojjTo. 


In  Eq.  3.70,  the  first  term  Rg  is 

Rg  = 

k=0 
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fc=0 
1  ^ 

=  ^M-j{mod{i  +  k-l,M  +  fc)(0d  -  <t>d)  +  INT{{i  -  1)/(M  +  k)){rd  -  0d)}] 

k=0 

=  R-s, 

because  (j)d  and  cancel  out  causing  the  exponential  terms  to  zero.  Thus  exponential  becomes 
equal  to  one  and  thus  the  sum  adds  to  K.  Also  R_,-  =  R^.  Therefore  Eq.  3.70  is 

Rg  =  Rg  +  Rj 
1  ^ 

+  Rgj  E  exp[-i{mod(i  +  k-l,M  +  k){<Pd  -  <l>j)  +  INT{{i  -  1)/(M  + 

k=l 
1  ^ 

+  Rjs;^  E  exp[-i{mod(i  +  k-l,M  +  k){4>j  -  <j>d)  +  INT{{i  -  1)/(M  +  k)){<^°j  - 

fc=i 

The  crosscorrelation  parts  (the  third  and  fourth  parts  of  above  equation)  go  to  zero  as  the  number  of 
smoothing  steps,  K,  goes  to  infinity.  Therefore  the  signals  are  progressively  decorrelated.  The  rate 
which  the  cross-correlation  part  approaches  zero  depends  on  the  hopping  frequency  and  propagation 
delay  (Reddy  et  al,  1987). 

The  smoothing  covariance  matrix  (Eq.  3.66)  is  used  to  calculate  the  weight  vector  Wq  in  the 
GSC,  as  Eq.  3.44 

Wa(n  d-l)  =  Wa(n)  +  ^cf  R(Wg  -  CoWa(Tl)).  (3.73) 

3.4-3  Summary.  A  jammer  that  is  correlated  with  the  desired  signal  arises  due  to  mul¬ 
tipath  propagation  or  intelligent  jamming  and  creates  a  special  problem  for  conventional  adaptive 
beamforming.  The  goal  of  this  section  was  to  illustrate  how  spatial  smoothing  combined  with  a 
LMS  adaptation  of  array  weights  can  remove  correlated  interference  in  a  GSC.  The  crosscorrelation 
between  desired  signal  and  jammer  signal  goes  to  zero  as  the  number  of  smoothing  steps,  K,  goes 
to  infinity. 
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3.5  Conclusion 


This  chapter  illustrated  the  theory  of  adaptive  signal  processing.  The  temporal  and  spatial 
filters  mathematical  formulations  and  procedures  for  their  solution  are  indeed  similar,  even  though 
they  are  used  in  different  application  areas.  Adaptive  temporal  filter  converges  to  the  optimum 
Wiener  solution  in  a  stationary  environment.  The  adaptive  beamformer  as  the  spatial  filter  forms 
a  weighted  combination  of  the  outputs  of  spatially  separated  antennas  in  order  to  suppress  propa¬ 
gating  signals  that  originate  from  different  locations.  The  LCMV-GSC  beamformer  is  used,  which 
minimize  the  output  power  with  related  to  the  undesired  signals.  The  adaptive  arrays  when  the 
interference(s)  is  correlated  with  the  desired  signal  were  totally  destroyed.  As  the  solution.  Section 
3.4  suggested  the  spatial  smoothing  technique  to  reduce  the  crosscorrelation  between  the  desired 
signal  and  the  jamming  signal.  The  next  chapter  will  simulate  the  jamming  suppression  in  the 
several  cases. 
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IV.  Simulation 


This  chapter  considers  the  simulation  of  jamming  suppression  in  several  cases.  A  LAS  antenna 
array  has  just  M  antenna  elements  [Fig.  1.1]  while  TDL  antenna  array  has  L  delayed  taps  in 
each  of  the  mth  antenna  element  channel  [Fig.  2.4].  Section  4.1  compares  the  performance  of 
wideband  interference  suppression  in  these  two  different  antenna  arrays.  Next  section  4.2  presents 
two  techniques  in  the  frequency  hopped  environment.  The  first  technique  is  the  conventional 
method  in  which  only  the  center  frequency  of  the  desired  signal  is  used  in  the  constraint  matrix 
(Eq.  3.19).  The  second  one  uses  each  hopping  frequency  of  the  desired  signal.  For  all  of  the  above 
cases,  we  assume  that  the  desired  signal  and  the  jamming  signals  are  uncorrelated.  In  section  4.3, 
we  investigate  the  case  when  the  desired  signal  and  jammer  signals  are  correlated.  Finally  the 
results  of  using  an  adaptive  filter  in  the  time  domain  for  narrowband  jamming  signals  suppression 
is  presented. 

4-1  Wideband  Jamming  Suppression  in  Antenna  Arrays, 

4.1.1  LAS  Antenna  Array  and  Frequency  Characteristics.  As  described  in  Section  3.3.2, 
the  signal  output  of  a  LAS  antenna  array  with  M  antennas  is 

yl^n)  =  W®^u(n).  (4.1) 

where  y{n)  is  the  output  of  antenna  array,  and  w  is  weight  vector  and  u(n)  is  input  signal  vector, 

W  =  [wi,W2,---  ,Wm] 

u{n)  =  [ui{n),U2{n),---UMin)] 

Ui{n)  =  u[nTs  —  mod{i  —  l,M)Tp\,  1  <i  <  M 
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where  Tp  =  ^  sin(^p)  is  the  interelement  propagation  delay  :  6p  is  the  direction  of  the  incident  signal. 
The  subscript  d  for  the  desired  signal  and  j  for  the  jammer  signal,  d  is  the  antenna  interelement 
space  which  is  a  half  wavelength  at  the  center  frequency.  Thus 


Tp 


—  sin  Oj, 
c 

—  sin  9„ 
2c 


The  signal  vector,  u(n)  can  be  written  a.s 

u(n)  =  Urf(n)  +  Uj  (n),  (4.2) 

where  Ud{n)  and  Uj(n)  are  vectors  containing  the  desired  and  interference  signals  respectively.  Ud 
and  uy(n)  are  assumed  to  be  uncorrelated. 

To  calculate  the  weight  coefficient  of  the  LAS  antenna  array,  a  linear  constraint  minimum  vari¬ 
ance  generalized  sidelobe  canceller  (LCMV-GSC)  beamformer  is  used  as  described  in  Section  3.3.5. 
The  desired  signal  direction  gain-only  constraint  equation  is  expressed  as, 

C^w  =  f,  (4.3) 

where  C  is  the  constraint  matrix,  C  =  d(0d,(^c}  and  f  is  the  response  vector,  f  =  1. 

d(6'd,Wc)  =  [1,  •••,  ^  (44) 

where  6d  is  the  incident  direction,  measured  with  respect  to  the  boresight  (i.e.,  the  normal  to  the 
array),  Wc  is  the  center  frequency  of  the  desired  signal. 

The  non-adaptive  beamformer  component  in  the  generalized  sidelobe  canceller  (GSC)  is  rep- 


4-2 


resented  by,  as  in  Eq.  3.36 


Wg  =  Cq=C(C"C)“^f.  (4.5) 

As  shown  in  Fig.  3.2,  the  weight  can  be  expressed  by  Eq.  3.37, 

w  -Wg-  CflWa,  (4.6) 

where  the  adaptive  beamformer  component  in  the  GSC  is  updated  as  in  Eq.  3.44, 

Wo(n  +1)  =  Wa(n)  +  /tx(n)e*(n) 

=  Wa(n)  +  /tCf  u(n)(wf  u(n)  -  wf  Cf  u(n))* 

=  Wa(n)  +  /tCf  u(n)u(n)^(wg  -  C„Wa(n))  (4.7) 

and  optimum  weight  vector  of  the  adaptive  beamformer  component  is  as  in  Eq.  3.48, 

Wa  =  (CfR„Ca)-^Cf  (4.8) 

The  output  frequency  transfer  function  is  (Wang  et  al,  1993) 

M-l 

H{(j}  =  (4.9) 

fc=0 

where  Tp  is  the  propagation  delay  time  between  the  elements  for  the  arriving  wave.  The  propagation 
delay  time  Tp  might  be  changed  to  Td  for  the  desired  signal  and  to  Tj  for  the  jammer  signal.  Hence 
the  output  frequency  transfer  function  of  the  LAS  antenna  array  for  the  desired  and  jamming 
signals  can  be  expressed  in  a  similar  form. 

If  the  direction  of  the  desired  signal  is  known,  the  weight  coefficients  Wk  for  the  antenna 
array  are  derived  by  means  of  the  LMS  adaptive  algorithm  as  in  Eq.  4.6  and  Eq.  4.7.  As  numerical 
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Frequency  Response  In  Wide  Jammer,  Jammer  BW  ■  4e-f04  Hz 


(a)  (b) 


Figure  4.1  (a)  Antenna  Array  Pattern,  (b)  Frequency  Characteristic  in  the  LAS  Antenna  Array 


examples,  the  antenna  pattern  and  the  frequency  characteristics  of  the  LAS  antenna  array  with  two 
elements  (M  =  2)  are  shown  in  Fig.  4.1  for  two  arriving  signals  (i  =  2).  In  the  antenna  pattern,  the 
desired  signal  incident  at  an  angle  of  30'’  from  boresight  is  passed  undistorted,  while  the  jamming 
signal  at  an  angle  of  -30*’  is  rejected.  This  is  more  clear  in  the  frequency  characteristics  graph  of 
a  Fig.  4.1.  The  solid  line  is  the  desired  signal  normalized  to  zero.  This  signal  passes  through  the 
antenna  array  whereas  the  dashed  line  shows  an  attenuation  of  the  jamming  signal.  This  figure 
also  shows  that  the  possible  blocking  bandwidths  of  the  jamming  signal  are  confined  to  narrow 
bandwidths. 


4.1.2  TDL  Antenna  Army  and  Frequency  Characteristics,  Fig.  2.4  shows  the  configuration 
of  the  TDL(tapped  delay  line)  antenna  array  in  which  several  tap  weights  are  provided  for  each 
antenna  channel.  In  adaptive  beamforming  for  narrowband  jammers,  the  directivity  is  constructed 
by  means  of  the  spatial  correlation  characteristics  of  the  signals  between  elements.  In  the  case 
of  the  TDL  array,  in  addition  to  the  capability  of  forming  an  identical  directivity,  the  temporal 
correlation  characteristics  of  the  signals  among  the  taps  can  be  used  in  such  a  way  as  to  improve 
the  frequency  characteristics  for  the  arriving  signal  with  a  broad  bandwidth. 
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The  signal  output  of  a  TDL  linear  array  with  M  antennas,  each  followed  by  a  tapped  delay-line 
with  L  taps,  is 

y{n)  =  w®u(n).  (4-10) 

where 


w  =  [wi,W2,---  ,wml] 

u(n)  =  [ui{n),U2in),---UML(n)] 


T  —  T 


Ui{n)  =  u[nTs  —  mod(i  —  l,M)Tp  —  INT{—j^)To] 


l<i<  ML 


and  where  Tp  is  the  propagation  delay  time  between  the  elements  for  the  arriving  wave.  The  delay 
time  To  of  each  tap  in  the  TDL  antenna  array  is  determined  by 


T - — 


(4.11) 


where  Rs  is  the  sampling  frequency  of  the  arriving  signal. 

In  the  constraint  equation  for  the  LCMV  beamformer  as  in  Eq.  3.19,  the  desired  signal 
direction  gain-only  constraint  is, 

f  =  [l,0...,0f  (4.12) 

^  d{0d,(^c)  Om  •••  Om 
Om  d(0d,Wc)  Om 


C  = 


(4.13) 


V  Om  Om  •••  d{0d,(Jc)/ 

where  d(0df^c)  is  M  X  1,  where  Oj  is  the  direction,  measured  with  respect  to  the  boresight,  Wc  is 
the  center  frequency  of  the  desired  signal,  and  is  expressed  as. 


d(0d,c^c)  =  ll, 


(4.14) 
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Antenna  response,  jammer  degree  -'30 


(a) 


Frequency  Response  in  Wide  Jammer,  Rj  -2e-i-04 


(b) 


Figure  4.2  (a)  Antenna  Array  Pattern,  (b)  Frequency  Characteristic  in  the  TDL  Antenna  Array 


and  the  column  vectors  Om  contain  M  zeros.  The  A-dimensional  gain  vector  f  describes  the 
frequency  response  of  the  beamformer  to  a  signal  impinging  on  the  array  from  the  desired  look 
direction  .  Thus  the  constraint  matrix  C  is  ML  x  L  dimensional  (Buckley  and  Griffiths,  1988). 

Using  the  weight  coefficient  which  is  updated  by  Eq.  4.6  and  Eq.  4.7,  the  output  frequency 
transfer  function  is  given  by 


M  L 

m=l  /=1 


(4.15) 


The  frequency  characteristics  of  the  TDL  antenna  array  output  with  two  elements  each  by 
having  five  delay  taps  are  shown  in  Fig.  4.2(b)  for  the  two  arriving  signals.  The  antenna  array 
pattern  of  the  computed  TDL  antenna  array  is  similar  to  those  in  Fig.  4.1(a).  The  output  frequency 
characteristics  pass  the  desired  signal  while  the  frequency  components  of  the  wideband  interference 
signal  within  the  bandwidth  can  be  sufficiently  suppressed.  As  a  result,  the  TDL  adaptive  antenna 
array  is  more  effective  for  the  suppression  of  wideband  interferences  than  the  LAS  antenna  array. 
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Item 

Desired  Signal 

Jamming  Signal 

Direction  at  Arrival 

30° 

-30° 

Carrier  Frequency 

15  X 

15  X  W{Hz) 

Bandwidth 

6  X  lO'^iHz) 

6  X  10“ 

Sampling  Frequency 

12  X  10“  Hz 

No.  of  TDL  array  elements 

8 

No.  of  TDL  tap  elements 

10,  20 

TDL  array  tap  delay 

1/Rs 

Table  4.1  Specification  for  First  Computer  Simulation 


4- 1-3  Comparison  of  Two  Antennas.  In  the  previous  section,  the  output  frequency 
characteristics  of  the  LAS  and  TDL  adaptive  antenna  array  were  studied.  In  this  section,  the 
performance  of  two  antenna  arrays  is  compared  using  output  signal  to  jammer  ratio  (SJR)  in 
the  presence  of  wideband  jammers.  The  wideband  signals  shown  in  Table  4.1  were  used  for  this 
simulation.  Their  fractional  bandwidth  is  a  j|o  =  0.04.  Fractional  bandwidth  larger  than  0.02  are 
defined  as  wideband  signals  [section  3.1.3]. 

The  performance  of  this  system  is  demonstrated  by  exploring  the  input  and  output  SJR  and 
improvement  defined  as, 


SJRo  = 

Jo 

Cl 

SJRi  = 

Si 

Ti 

(4.16) 

SJRo 

IMP  = 

SJRi 

(4.17) 

where  Si  and  So  is  the  desired  signal  input  and  output  power,  and  Ji  and  Jo  is  the  jamming  signal 
power,  respectively. 

The  LMS  weight  estimate  is  updated  at  each  time  n  by  Eq.  4.6  and  Eq.  4.7.  Fig.  4.3  shows 
the  IMP  of  the  LAS  and  TDL  antenna  array  a.s  the  function  of  iteration.  Both  systems  have  a 
good  convergence  characteristic.  IMP  of  the  TDL  antenna  array  converges  to  higher  value  than 
IMP  of  the  LAS  antenna  array.  In  order  to  see  this  more  clearly,  the  IMP  is  separated  into  the 
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Improvement  as  the  Function  of  Iteration 


Figure  4.3  IMP  in  LAS  and  TDL  Antenna  Array  using  LMS  Weight  Estimate 

desired  signal  power  {So  and  Si)  and  the  jamming  signal  power  {Jo  and  Ji).  Fig.  4.4(a)  shows  the 
desired  signal  power  variation  as  a  function  of  iteration.  The  solid  line  is  the  desired  signal  input 
power  {Si),  and  the  dashed  line  is  the  desired  signal  output  power  {So)  of  the  LAS  antenna  array. 
The  dotted  line  is  So  of  the  TDL  antenna  array.  The  desired  signal  output  power  (S'o)  is  close  to 
the  signal  input  power  {Si)  in  both  antennas  while  the  jammer  signal  is  suppressed  significantly 
[Fig.  4.4(b)].  In  Fig.  4.4(b),  the  dotted  line  shows  the  jamming  signal  output  power  in  the  LAS 
antenna  array  while  the  dashed  line  shows  it  in  the  TDL  antenna  array.  As  a  result,  Fig.  4.4(b) 
shows  that  the  jamming  signal  is  suppressed  more  efficiently  in  the  TDL  antenna  array  than  in 
the  LAS  antenna  array  as  much  as  lOdB.  Therefore,  we  can  prove  the  TDL  antenna  array  is  more 
efficient  than  the  LAS  antenna  array  for  wideband  jammers  suppression.  Fig.  4.4(c)  and  Fig.  4.4(d) 
show  the  same  effect  as  Fig  4.1(b)  and  Fig  4.2(b). 

The  next  experiment  shows  the  effect  of  increasing  the  delayed  taps  from  10  to  20  for  a  fixed 
number  of  antenna  elements  (M=8).  As  expected,  the  performance  of  the  TDL  antenna  array  is 
improved  by  increasing  the  delay  taps  from  10  to  20  [Fig.  4.5].  The  desired  signal  power  (5o)  is  not 
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Desired  Signal  Power  (dB) 


Signs!  Power  (d^ 


Item 

Desired  Signal 

Jamming  Signal 

Direction  at  Arrival 

SO" 

-30° 

Carrier  Frequency 

15  X  lO^iHz) 

15  X  10^  (Hz) 

Bandwidth 

3  X  10^ (Hz) 

Sampling  Frequency 

12  X  10*^  Hz 

Table  4.2  Specification  for  Second  Computer  Simulation 


changed  much  when  compared  to  the  input  signal  power  (Si)  [Fig.  4.6(a)].  On  the  other  hand,  as 
the  number  of  delay  taps  increases,  the  jamming  signal  power  is  reduced  significantly  [Fig.  4.6(b)]. 

The  next  simulation  considers  the  changes  when  the  fraction  bandwidth  increased  form  0.04 
to  0.2  as  shown  in  Table  4.2.  This  means  that  the  BW  of  signals  in  specification  of  computer 
simulation  is  changed  from  6  x  10^  to  3  x  10^.  Fig.  4.7  still  shows  the  TDL  antenna  array  is  more 
effective  in  suppression  of  the  jamming  signal  than  the  LAS  antenna  array.  If  we  see  Fig.  4.8(a), 
however,  there  are  difference  from  Fig.  4.6(a).  The  desired  signal  output  power  (So)  is  decreased  in 
the  both  antenna  arrays.  This  is  due  to  the  fact  that  the  LCMV-GSC  beamformer  constrained  the 
direction  and  center  frequency  of  the  desired  signal  using  a  constraint  matrix.  If  the  input  signal 
is  narrowband,  there  is  no  distortion  of  the  desired  signal.  In  this  experiment,  however,  the  input 
signal  has  extremely  wide  bandwidth  (fractional  BW  =  0.2).  As  the  signal  bandwidth  increases, 
the  distortions  are  increased  in  the  desired  signal  except  the  center  frequency  used  in  the  constraint 
matrix.  As  you  see  Fig.  4.4  (c)  and  (d),  the  desired  signal  has  unity  gain  at  the  center  frequency. 
Other  frequencies  are  attenuated.  Thus,  if  the  desired  signal  has  wide  bandwidth,  this  signal  is 
distorted  passing  the  antenna  arrays. 

Another  interesting  point  in  Fig.  4.8(a)  is  that  the  desired  signal  output  power  (S'o)  of  the 
LAS  antenna  array  is  slightly  higher  than  the  So  of  the  TDL  antenna  array.  The  phase  difference 
between  the  weight  coefficients  depends  on  the  time  delay  and  frequencies  of  the  desired  signal. 
The  constraint  however  fixes  the  phase  with  respect  to  the  center  frequency  and  the  time  delay. 
Due  to  difference  in  phase  between  the  real  weights  and  the  constraint,  distortion  occurs.  In  the 
LAS  antenna  array,  only  M  antenna  elements  affect  the  distortion  while  the  distortion  is  affected 
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Improvement  as  the  Function  of  Iteration 


Figure  4.7  IMP  in  LAS  and  TDL  Antenna  Array  using  LMS  Weight  Estimate 

by  not  only  M  antenna  elements  but  also  (M  —  1)L  delayed  taps  in  the  TDL  antenna  array. 

Fig.  4.8(b)  indicates  the  jamming  signal  input  and  output  powers  as  a  function  of  iteration. 
Both  systems  show  the  good  convergence  for  jamming  signal  suppression.  As  the  number  of  the 
iteration  increases,  the  TDL  antenna  array  suppresses  the  jamming  signal  more  than  the  LAS 
antenna  array. 

The  next  experiment  shows  the  performance  effect  by  increasing  the  delayed  taps  from  10  to 
20  in  the  extreme  wideband  signals.  The  IMP  of  the  TDL  antenna  array  decreases  even  though 
the  number  of  delayed  taps  increases,  as  shown  in  Fig.  4.9.  The  reason  is  that  the  desired  signal 
distortion  was  more  serious  with  the  larger  number  of  delayed  taps  [Fig  4.10(a)].  Furthermore  the 
capability  of  the  jamming  signal  suppression  is  not  even  changed  [Fig  4.10(b)]. 
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signal  Power  (dB) 


Improvement  as  a  Function  of  Iteration 


Figure  4.9  The  Comparison  Between  K  =  10  and  iT  =  20  in  the  TDL  Antenna  Array 

Input  and  Output  Deeired  Signal  Power  in  Iteration  Input  and  Output  Jammer  Signal  Power  In  Iteration 


(a) 


(b) 


Figure  4.10  The  Comparison  Between  K  =  10  and  if  =  20  in  the  TDL  Antenna  Array  (a)  The 
Desired  Signal  Part,  (b)  The  Jamming  Signal  Part 
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Jamming  Suppression  in  Frequency- Hopped  Environment 


This  section  examines  the  output  SJR  performance  of  adaptive  arrays  in  a  Frequency-Hopped 
(FH)  environment.  The  combined  effects  of  the  arrays  interference  nulling  ability  and  the  interfer¬ 
ence  protection  of  FH  signals  result  in  a  more  robust  interference  rejection  communication  system. 
As  concluded  in  the  previous  section,  the  TDL  antenna  array  should  be  more  effective  than  the 
LAS  antenna  array  for  frequency-hopped  jamming  signal  suppression,  because  FH  signals  also  have 
wide  bandwidths. 

In  the  next  example,  two  different  techniques  are  used  for  frequency  hopped  jamming  sup¬ 
pression  in  the  TDL  antenna  array.  The  first  conventional  technique  is  the  same  as  the  method  of 
section  4.1.2,  which  uses  the  center  frequency  in  the  constraint  matrix.  The  second  method  uses  a 
hopping  frequency  in  the  constraint  matrix. 

4-2.1  Conventional  Technique  for  Frequency-Hopped  Environment.  For  the  experiment, 
the  TDL(tapped  delay  line)  antenna  array  is  used.  With  M  antenna  elements  spaced  a  half  wave¬ 
length  apart  at  the  center  frequency  of  the  FH  band,  Wc  =  ^irfc.  A  desired  signal  is  incident  on 
the  array  from  direction  6^  and  with  hopping  frequencies  u)h  =  27r/*,  and  an  jamming  signal  from 
direction  6j  and  at  its  hopping  frequencies  tOj  =  2nfj.  The  input  signal  vector,  u(ti)  can  be  written 
as, 

u(n)  =  Ud(n)  -t-  Uj(n),  (4-18) 

where  Ud(n)  and  Uj(7i)  are  vectors  containing  the  desired  and  interference  signals  respectively.  Ud 
and  Uj(n,)  are  assumed  to  be  uncorrelated. 

The  delay  time  To  of  each  tap  in  the  TDL  antenna  array  is  determined  by 


To 


(4.19) 
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where  fc  is  the  carrier  frequency  of  the  arriving  signal.  R.T  Compton  shows  the  coefficient  7  in  his 
journal  (Compton,  1988)  as 

;^  <  7  <  |,  (4.20) 

where  B  is  determined  as  the  normalized  jammer  bandwidth, 

Jc 

The  form  of  the  desired  signal  vector  and  interference  signal  vector  are 

u(n)  =  u{n)e{h),  Th-i  <n<Th,  (4-21) 

where  Th-i  <  n  <Th  is  the  time  duration  of  the  hop  interval,  and  the  signal  arrival  vector  is, 

ei{h)  =  exp[-i{mod(i  -  l,M)^p(h)  -  7iVr(^)<A;(h)}].  (4.22) 

Here,  <l)d{h)  is  the  desired  signal  interelement  phase  shift  during  hop  h  and  (l>j  is  the  interference 
signal  interelement  phase  shift  during  hop  j,  i.e, 

<Pd{h)  =  UhTd 
<l>j{h)  = 

and  <j>p{h)  is  the  phase  shift  due  to  taps  delay  time  To,  so, 

rM  =  cvhTo 
r^ih)  =  cdjTo. 

The  constraint  matrix  and  the  weight  coefficient  of  the  TDL  antenna  array  is  the  same  as  the 
section  4.1.2.  The  frequency  characteristics  of  the  frequency  hopped  signals  and  the  comparison 
between  the  LAS  antenna  array  and  the  TDL  antenna  array  are  very  similar  to  Fig.  4.2  and  Fig.  4.3, 
respectively  since  frequency  hopped  signals  belong  to  the  wideband  signals. 
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4.2.2  New  Technique  for  Frequency- Hopped  Environment.  The  signal  output  of  the  TDL 
linear  array  with  M  antennas,  each  followed  by  a  tapped  delay-line  with  L  taps,  is 


where 


y(n)  =  w"u, 


(4.23) 


W  =  [WI,W2,---  ,Wml] 

u  =  u(n)e{h),  Th-\  <n<Th. 

The  signal  arrival  vector  is 


ei{h)  =  exp[-j{mod{i  -  l,M)^p(h)  -  7iVr(^)^;(h)}].  (4.24) 

Here,  <j}d{h)  is  the  desired  signal  interelement  phase  shift  during  hop  h,  and  <j>i  is  the  interference 
signal  interelement  phase  shift.  And  4>p{h)  is  the  phase  shift  due  to  delay  time  To- 

In  the  constraint  equation,  the  desired  signal  direction  gain-only  constraint  is, 

f=[l,0...,0f,  (4.25) 


and 


/d{ed,o;h) 

Om 

Ok 

d(^'d,w/,)  . 

o  o 

\  Om 

Om 

■  d{ed,u}h)) 

(4.26) 
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Improvement  as  a  Function  of  Iteration 


Figure  4.11  Comparison  of  IMP  Between  Two  Different  Techniques 


where  d{6d,‘jJh)  is  Af  x  1,  where  9d  is  the  direction,  measured  with  respect  to  the  boresight,  Wh  is 
the  hopping  frequency  of  the  desired  signal  respectively,  and  is  expressed  as 


d(0d,wh)  =  [l, 


(4.27) 


The  difference  here  compared  to  the  conventional  method  is  that  the  hopping  frequency  is  used 
in  the  constraint  matrix  instead  of  the  center  frequency.  The  output  frequency  transfer  function 
remains  the  same  as  Eq.  4.15. 


4-2.3  Comparison  of  the  Two  Techniques.  The  frequency-hopped  signals  generated  by 
Table  4.2  are  used  for  this  experiment.  They  are  incident  into  the  antenna,  each  having  their  own 
direction.  The  TDL  antenna  array  weights  are  updated  by  the  LCMV-GSC  algorithm,  as  described 
in  Eq.  4.6  and  Eq.  4.7.  Fig.  4.11  shows  the  IMP  of  the  two  techniques  for  the  TDL  antenna  array. 
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As  a  result,  IMP  of  the  new  technique  is  higher  than  IMP  of  the  conventional  technique  as  much 
as  15dB. 

This  can  be  explained  in  Fig.  4.12.  In  Fig.  4.12(a),  the  desired  signal  output  power  of  the 
new  technique  is  same  as  the  input  power  while  So  of  the  conventional  method  is  reduced.  In  the 
constraint  matrix,  the  new  technique  constrains  each  hopping  frequency,  therefore  the  frequency 
characteristics  of  the  desired  signal  attempts  to  match  the  hopping  frequency.  When  the  hopping 
frequency  of  the  desired  signal  is  changed,  the  frequency  characteristic  of  the  desired  signal  has  a 
unity  gain  at  the  hopping  frequency  as  shown  in  Fig.  4.12(c).  Therefore,  there  is  no  distortion  in  the 
desired  signal  part.  On  the  other  hand,  in  the  conventional  method,  the  frequency  characteristics  of 
desired  signal  maintains  the  center  frequency  with  a  unity  gain  even  though  the  signal  frequencies 
are  changed  to  the  other  hopping  frequencies.  Thus  there  is  some  distortion  at  the  other  hopping 
frequencies  [Fig.  4.12(d)].  Fig.  4.12(b)  shows  that  the  jamming  suppression  capability  of  the  new 
technique  is  greater  than  the  conventional  technique. 

In  the  new  technique,  the  next  experiment  investigated  the  change  when  the  the  number  of 
taps  increases  from  10  to  20.  Fig  4.13  shows  the  IMP  for  two  different  numbers  of  delay  taps. 
Unlike  the  conventional  method,  the  IMP  increases  by  increasing  the  number  of  taps.  This  occurs 
because  the  desired  signal  power  (So)  is  not  attenuated  while  the  jamming  signal  is  more  suppressed 
by  increasing  the  number  of  delay  taps  (Fig.  4.14(a)(b)). 
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Deeired  Signal  Power  (dB) 


Input  and  Output  Desired  Signal  Power  In  Iteration  input  and  Output  Jammer  Signal  Power  as  a  Function  of  Iteration 


Iteration  Iteration 


Figure  4.12  Input  and  Output  Power  of  New  and  Conventional  Technique  (a)  Desired  Signal  Part 
(b)  Jamming  Signal  Part,  Frequency  Characteristics  of  (c)  New  Technique  and  (d) 
Conventional  Technique  Respectively 
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Desired  Signal  Power  (dB) 


Improvement  as  Function  of  Iteration 


Figure  4.13  The  Comparison  Between  iT  =  10  and  IT  =  20  in  the  New  Technique 


Input  and  Output  Desired  Signal  Power  In  Iteration  Output  Jamming  Signal  Power  as  a  Function  of  Iteration 


(a) 


(b) 


Figure  4.14  The  Comparison  Between  FT  =  10  and  FT  =  20  in  the  New  Technique  (a)  The  Desired 
Signal  Part,  (b)  the  Jamming  Signal  Part 
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4-3  The  Performance  of  a  Spatial  Smoothing  Technique  for  Correlated  Signals 


This  section  presents  an  example  case  in  which  the  desired  signal  and  the  jamming  signal  are 
correlated.  The  simulation  uses  the  GSC  adaptive  weights,  using  the  method  which  is  described 
in  section  3.4.  Both  signals  have  wide  bandwidths,  and  the  correlation  coefBcient  between  desired 
signal  and  jammer  signal  is  1, 


P  = 


(4.28) 


where  Raj  is  the  cross  correlation  between  the  desired  signal  and  the  jamming  signal  and  Rs,Rj 
are  the  autocorrelation  of  the  desired  signal,  the  jamming  signal  respectively. 

For  this  experiment,  the  LCMV  beamformer  with  the  TDL  antenna  array  is  used  as  in  the 
previous  section.  The  input  signals  and  the  weight  vectors  are  the  same  as  in  Section  4.1.  The  only 
difference  from  section  4.1  is  that  the  signals  are  correlated  with  each  other. 

The  correlation  matrix  used  for  updating  the  weights  is,  as  in  Eq.  3.66, 


1  ^ 


fc=i 


(4.29) 


The  sub-array  size  M  was  fixed  at  4  and  the  number  of  delayed  taps  is  10.  The  number 
of  smoothing  steps  K  is  5.  The  desired  signal  input  power  and  the  jamming  signal  input  power 
are  equal.  Fig.  4.15  shows  the  SJR  of  three  different  cases  as  a  function  of  iterations  .  The 
solid  line  shows  the  SJRg  when  the  signals  are  uncorrelated,  while  the  dashed  line  and  the  dot- 
and-dashed  line  show  the  correlated  signals  case.  The  dashed  line  displays  before  decorrelation, 
which  means  K  =  1.  At  that  time,  the  SJRq  is  reduced  significantly,  because  the  antenna  can’t 
put  a  null  to  the  direction  of  the  jamming  signal  as  shown  in  Fig.  4.16(c).  The  jamming  signal 
is  a  replica  of  the  desired  signal.  The  antenna  array  accepts  the  jamming  signal  as  the  desired 
signal.  This  fact  shows  that  the  adaptive  antenna  array  is  useless  in  the  case  of  smart  jammers 
and  multipath  signals  because  those  jamming  signals  and  multipath  signals  are  highly  correlated 
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Output  SJR  in  iteration 


Figure  4.15  The  Comparison  of  SJR  After  Performing  the  Spatial  Smoothing  Technique  in  the 
Presence  of  Correlated  Signals 

with  the  desired  signal.  However,  after  performing  the  smoothing  technique  to  reduce  the  cross¬ 
correlation  between  the  desired  signal  and  jammer  signal,  the  performance  of  adaptive  antenna 
array  improved  significantly  as  shown  in  Fig.  4.15  by  the  lines  as  much  as  20dB. 

The  desired  signal  power  variation  is  shown  in  Fig.  4.16(a)  as  a  function  of  iteration.  As 
explained  earlier,  the  desired  signal  has  some  distortion.  The  output  power  of  the  desired  signal 
before  decorrelation  has  a  higher  value  than  even  the  uncorrelated  case.  This  occurs  because  some 
of  the  jamming  signal  is  added  to  the  desired  signal.  The  jamming  signal  output  power  (Jo) 
before  decorrelation  signal  is  close  to  the  jammer  input  power  ( Jj)  (Fig.  4.16(b)),  which  shows  no 
suppression  at  all  in  the  jamming  signal  (Fig.  4.16(c)).  After  performing  the  smoothing  technique 
however,  the  jamming  signal  output  power  is  reduced,  which  shows  the  smoothing  technique  works 
for  suppression  of  the  jamming  signal  in  the  adaptive  antenna  array.  Fig  4.16(d)  shows  the  frequency 
characteristics  after  performing  the  spatial  smoothing  technique. 
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Signal  Power  (dB) 


Input  and  Output  Signal  Power  In  Iteration 


Iteration 


Input  and  Output  Jammer  Signal  Power  In  iteration 


Iteration 


Figure  4.16  (a)  The  Desired  Signal  Power  Variation,  (b)  The  Jamming  Signal  Power  Variation, 

(c)  and  (d)  Frequency  Characteristics  Before  Decorrelation  and  After  Decorrelation, 
Respectively 
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4-4  Adaptive  Filter  for  Narrowband  Interference  Suppression 


The  desired  signal  used  here  is  the  same  as  that  which  is  generated  in  Section  4.1.  The  desired 
signal  is  assumed  to  be  known  in  the  receiver.  The  jamming  signal  is  a  single  tone  jammer.  The 
number  of  taps  is  12.  The  input  signal  vector,  u  and  weight  vector  w  are, 

u  =  [u{n),u{n—\) . M(n  —  M  +  1)]^ 

w  =  [W1,W2 . wm]^, 

where  M  =  12  is  the  number  of  taps. 

The  instantaneous  estimated  value  of  the  tap-weight  vector  at  time  n  -f  1  is  computed  by 
using  the  simple  recursive  relation  as  in  Eq.  3.11, 

w(n  +  1)  =  w(n)  -t-  /fu(n)[ci*(n)  -  u'^(n)w].  (4.30) 

where  d{n)  is  the  desired  signal,  and  p  =  is  the  step  size  parameter.  \max  is  the  largest 

eigenvalue  of  correlation  matrix  R. 

The  performance  of  this  adaptive  filter  is  demonstrated  by  exploring  the  SJR  improvement. 
Fig.  4.17(a)  shows  the  SJR  improvement  as  afunction  of  iterations.  Fig.  4.17(b)  shows  the  frequency 
characteristics  of  the  adaptive  filter.  This  system  inserts  the  notch  filter  at  the  jammer  carrier 
frequency  at  the  cost  suppressing  that  frequency  in  the  desired  signal.  Therefore  this  system  only 
works  for  narrowband  jamming  signal  suppression.  The  jamming  signal  is  suppressed  well  after 
just  a  few  first  iterations  [Fig.  4.17(d)].  On  the  other  hand,  the  desired  signal  is  also  suppressed  at 
the  first  iterations.  However  as  the  number  of  iterations  increases,  the  desired  signal  recovers  its 
input  power  [Fig.  4.17(c)].  This  means  the  notch  filter  is  narrowing  in  order  to  avoid  suppressing 
the  desired  signal.  Eventually,  the  desired  signal  output  power  approaches  the  input  power.  Thus 
the  SJR  improvement  of  this  system  reaches  the  amount  of  jamming  signal  suppression. 
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Adaptive  Frequency  Response 


The  number  of  Tap.  -  12 
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Figure  4.17  (a)  Improvement  of  Adaptive  Filter,  (b)  Frequency  Characteristics  of  Adaptive  Filter, 

(c)  Si  and  Sg,  (d)  Ji  and  Jg 
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4-5  Conclusion 


This  chapter  considered  the  simulation  of  jamming  signal  suppression.  The  TDL  antenna 
array  is  more  effective  than  the  LAS  antenna  array  to  suppress  the  wideband  interference  in  the 
LCMV-GSC  beamformer.  In  the  frequency-hopped  environment,  a  new  technique,  which  the  hop¬ 
ping  frequencies  are  used  in  the  desired  signal  direction  gain-only  constraint  matrix,  suppressed 
wideband  jamming  signal  more  efficiently  than  the  conventional  method.  The  smoothing  technique 
was  demonstrated  to  reduce  the  cross-correlation  of  the  desired  and  jamming  signal.  After  perform¬ 
ing  this  technique,  the  jamming  signal  was  suppressed  significantly.  Finally,  the  performance  of 
the  adaptive  filter  was  shown  for  suppression  of  the  narrowband  jamming  signal.  The  next  chapter 
will  conclude  this  research  and  recommend  research  topics  for  future  study. 
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V.  Conclusion  and  Recommendations  for  Future  Research 

5.1  Conclusion 

Interference(s)  is  a  significant  problem  in  military  communication  operations.  This  thesis 
examined  the  suppression  of  interference  using  the  SJR  as  a  metric  for  comparing  the  adaptive 
temporal  filter  and  adaptive  spatial  filter  (i.e.  adaptive  beamforming).  In  these  systems,  the  opti¬ 
mum  weights  are  selected  to  optimize  the  response  based  on  the  statistics  of  the  data.  The  data 
statistics  are  often  unknown  and  may  change  with  time  so  adaptive  algorithms  are  used  to  obtain 
weights  that  converge  to  the  statistically  optimum  solution.  Here,  the  most  well-known  and  simple 
algorithm,  LMS  algorithm  was  used. 

For  the  class  of  beamformer,  a  procedure  was  presented  for  the  design  and  implementation  of 
a  LCMV  beamformer  for  wideband  interference  suppression.  This  beamformer  attempts  to  achieve 
the  minimum  output  variance  subject  to  the  constraint  on  the  gain  in  the  desired  signal  direction. 
The  structure  used  extends  the  GSC  to  implement  and  to  analyze  the  LCMV  beamformer.  The 
CSC  separates  the  weight  vector  with  constrained  and  unconstrained  components.  The  constrained 
components  represent  the  gain  and  phase  of  the  desired  signal  and  its  direction.  The  unconstrained 
components  are  adjusted  using  the  LMS  adaptive  algorithm.  Using  the  LCMV-GSC  beamformer, 
it  was  shown  that  the  TDL  adaptive  antenna  array  is  more  effective  for  the  suppression  of  wideband 
jammer  suppression  than  the  LAS  adaptive  antenna  array. 

Further,  in  the  frequency-hopped  environment,  the  new  technique  was  introduced  which  used 
each  hopping  frequency  for  the  desired  signal  direction  gain-only  constraint  in  the  LCMV  beam- 
former.  This  new  technique  improved  the  capability  of  wideband  interference  suppression  over  the 
conventional  technique  which  used  only  the  center  frequency  in  the  constraint  matrix.  The  reason 
is  that  the  new  technique  does  not  distort  the  desired  signal,  while  the  conventional  technique 
severely  distorts  the  desired  signal. 

Then  the  effect  of  the  spatial  smoothing  in  the  presence  of  the  correlated  signals  was  ana- 
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lyzed.  When  the  desired  signal  was  correlated  with  the  jamming  signal,  the  LCMV  beamforming 
performed  poorly  because  the  adaptive  antenna  array  accepted  the  jamming  signal  as  the  desired 
signal.  Therefore  it  couldn’t  put  a  null  in  the  direction  of  the  jamming  signal.  After  performing 
the  spatial  smoothing  technique,  however,  the  adaptive  antenna  array  inserted  a  deep  null  for  sup¬ 
pression  of  the  jamming  signal. 

Finally,  the  adaptive  temporal  filter  was  shown  to  adequately  suppress  narrowband  jammers. 
The  drawback  of  this  system  however,  is  that  it  can  not  be  used  to  suppress  wideband  interference. 

5.2  Recommendation  for  Future  Study 

There  are  three  main  areas  of  research  that  can  be  derived  from  this  thesis  effort.  These  three 
areas  are, 

•  Adaptive  Beamforming  for  Wideband  Interference. 

•  Spatial  Smoothing  in  the  Presence  of  Correlated  Signals. 

•  Adaptive  Filter  for  Narrowband  Interference. 

In  the  first  of  these  areas,  there  are  the  desired  signal  distortions  when  using  the  conventional 
technique.  Future  study  is  recommended  to  investigate  a  way  to  decrease  this  distortion.  Since  the 
desired  signal  is  distorted,  the  SJRo  decreases  even  though  the  number  of  delayed  taps  increases. 
Thus  it  is  recommended  to  study  the  optimum  number  of  delayed  taps  and  antenna  elements.  This 
thesis  was  limited  to  a  single  jammer  and  jamming  signal  direction.  The  performance  of  interfer¬ 
ence  suppression  needs  to  be  observed  for  a  larger  number  of  jammers  and  for  changing  directions. 
Because  of  the  computational  consideration,  the  partially  adaptive  beamformer  with  array  com¬ 
posed  of  large  numbers  of  sensors  has  been  investigated  and  remained  as  a  good  topic. 

Another  area  for  possible  research  is  suppression  of  the  jamming  signal  correlated  with  the 
desired  signal.  This  area  arises  due  to  multipath  propagation  or  smart  jamming  and  creates  a 
serious  problem  for  conventional  adaptive  beamformer.  The  spatial  smoothing  technique  is  based 
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on  averaging  to  destroy  the  correlation  between  the  signal  and  interference.  However,  correlated 
interference  suppression  still  remains  an  important  research  topic. 

A  final  recommended  area  for  future  research  regards  temporal  adaptive  filters.  It  is  recom¬ 
mended  to  compare  the  performance  of  interference  suppression  of  the  time  domain  adaptive  filter 
with  the  transform  domain  adaptive  filter  and  the  transform  domain  excisor. 
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Appendix  A.  Constrained  Optimization 

The  typical  constraints  optimization  problem  has  the  form 

min/(x),  subject  to  g(x)  =  0,  (A.l) 

where  f(-)  is  the  scalar-valued  objective  function  and  g(-)  is  the  vector-valued  constaint  function. 
The  classical  approach  to  solving  a  constrained  optimization  problem  is  the  method  of  Lagrange 
multipliers,  A.  The  Lagrangian  of  a  constrained  optimization  problem  is  defined  to  be  the  scalar¬ 
valued  function 

-^^(x,A)  = /(x)-|- Ag(x).  (A.2) 

In  the  case  at  hand  one  minimize  w^RiW  subject  to  the  linear  constraint  C-^w  =  f .  w  is  weight 
vector  and  Rj:  is  autocorrelation  vector  of  input  vector,  x.  The  constraint  g{x)  =  C^w  -  f  is  a 
scalar-valued  one  where  C  is  the  constraint  matrix  and  f  is  response  vector.  The  Lagrangian  is 

I,(w,  A)  =  ■w^Ra;W -h  A(C'^w  -  f).  (A.3) 

Its  gradient  is  2Rw  +  AC  with  a  solution, 


^  opt 


Aopf  R 


-1 


c 

2' 


(A.4) 


To  find  the  value  of  the  Lagrange  multiplier,  this  solution  must  satisfy  the  constraint.  Imposing 
the  constraint, 

AoptC'^'R-^C  =  -2f.  (A.5) 


By  solving  with  Xopt, 


Xopt  — 


2f 

C^R-iC 


(A.6) 
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The  total  solution  is 


Wopt 


fR-^C  ,  , 

CJ?R-1C 
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Appendix  B.  Matlab  Coding 


B.l  Generation  of  the  Signals 

Rs  =  $2*10*5$;  '/,  The  sampling  frequency. 

fl  =  50000;  y,  The  carrier  frequency  of  the  desired  signal. 

f2  =  50000;  y,  The  carrier  frequency  of  the  jamming  signal. 

N  =  5000;  y,  The  number  of  samples. 

B.1.1  General  Wideband  Signals. 

x=randn(l,N) ;  */,  Generate  the  random  noise 

fcutll=(Rc)/Rs*2;  '/,  Cut  frequency  of  the  filter. 

[n,wn]=cheb2ord(f  cutll  ,fcutll+0 .02,0.05,60)  ;  */,  Decide  order  of  the  filter. 

[B,A]=cheby2(n,60,wn) ;  '/,  Decide  the  filter  coefficient. 

signal=f  ilter(B,A,x) ;  "/,  Filters  the  data  by  B  and  A. 

for  i=l:N 

expon(i,l)=cos(2*pi*fl/Rs*(i-l));  ’/,  Carrier  frequency 

end 

signal=signal*expon; 

B.l. 2  Frequency  Hopping  Signals. 


y.  Generate  R  bits  bipolar  data  sequence  having  P  data  rate, 
bi  =  sign(randn(l,R)) ; 
data  =  bi’*ones(l ,P) ; 
datasequence  =  reshape (data’ , 1 ,P*R) ; 


7t  Generate  BFSK  signal  of  R  bits. 
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bfsk=zeros(M,length(n)) ; 


7,  Initial  value  of  BFSK  signal. 


wl=1000; 

w2=2000; 

n=l:20; 


7,  The  first  frequency  of  BFSK  signal. 

7,  The  second  frequency  of  BFSK  signal. 


bfskl=cos(2=t=pi*wl.*n./fsim) ; 
bfsk2=cos(2*pi*w2/fsiin.*n) ; 
for  i=l:R 
if  bi(i)==l 

bfsk(i , ; )=bf ski ; 
else 

bfsk(i, :)=bfsk2; 


7o  The  data  sequence  of  the  first  BFSK  signal. 

7«  The  data  sequence  of  the  second  BFSK  signal. 

70  Decide  the  bipolar  data  sequence 

71  using  BFSK  sequences. 


end 

end 

databfsk=reshape(bfsk’ ,l,M*length(n));  7»  Generate  the  BFSK  signal. 

7(  Generate  hopping  signal. 

h=zeros(pg,N) ;  %  Initial  value  of  hopping  signal, 

for  i=l:pg 

nuin(i)=round(rand(l , l)*pg) ; 
end 

for  i=l:pg; 

ffl(i)=fh*num(i)+fi;  7.  Decide  the  hopping  frequency, 

for  n=l:N; 

h(i,n)=cos(2*pi*ffl(i)*n/fsiin) ;  %  Make  the  hopping  sequence, 

end 


end 

h=reshape(h’ ,l,N*pg) ; 
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’/,  Generate  the  FH  (Frequency  Hopping)  signal. 


s ignal=dat  abf  sk . *h ; 


B.2  LCMV  Adaptive  Beamforming 

M  =  12; 

K  =  5; 

mul  =  0.0002; 
thd  =  30*pi/180; 
thi  =  -30*pi/180; 

Td  =  2/f l+sin(thd) ; 

Ti  =  2/f l*sin(thi) ; 
ffl 


y.  The  number  of  antenna  elements. 
*/,  The  number  of  delayed  Taps. 

%  The  step  size  parameter. 

7,  The  desired  signal  direction, 
y.  The  jamming  signal  direction. 

7,  Propagation  delay  by  signal. 

%  Propagation  delay  by  jammer. 

%  The  signal  hopping  frequency. 


y.  Generate  the  constraint  matrix  for  the  LMS  adaptive  array  antenna, 
for  ii=l:M 

d(ii, l)=exp(-sqrt(-l)*2*pi*(fc)/2/ (fc)*sin(thd)*(ii-l)) ; 
end 

C=d;  y,  Constraint  matrix  for  desired  direction  gain-only. 

g=l;  y.  The  response  of  the  beamformer. 

Wq=C*inv(C’*C)*g;  7,  Wq  is  the  quiescent  weight  vector. 


y,  The  blocking  Matrix. 

[m,n]  =  size(C); 

C  =  zeros(m,m) ; 

CQ,R]  =  qr(C); 

Ca  =  Q(:,n+l:m);  '/,  The  blocking  matrix. 
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°/o  Generate  the  constraint  matrix  for  the  TDL  adaptive  array  antenna, 
for  j=l:K 
for  ii=l:M 

d(ii,l)=exp(-sqrt(-l)*2*pi*(fc)/2/(fc)*sin(thd)*(ii-l)) ;  '/,  Array  response 
end  ’/,  vector. 

C( : ,  j)=[zeros(M*(j-l) ,  1)  ;  d;  zeros(M*(K-j)  ,1)]  ;  '/.Constraint 

end  '/,  matrix 

f=[l  zeros(l,K-l)]  ’ ;  '/,  The  response  of  the  beamformer. 

Wq=C*inv(C’*C)*g;  %  Wq  is  the  quiescent  weight  vector. 

[Ca]=block(C)  ;  '/,  Blocking  Matrix 

'/,  The  LMS  antenna  input  signal  vector  and  LMS  algorithm. 
si=0; 

for  i=Nl:N 
si=si+l 
for  j=l:M 

ud(j ,si)=inputd(i-Td*Rs*(j-l)) ; 
ui(j ,si)=winputi(i-Ti*Rs*(j-l)) ; 
end 

u(  :  ,si)  =  (ud( : , si)+ui( :  ,si) )  ;  ’/,  The  input  signal  vector. 

R=u( :  ,si)*u( :  ,si)  ’ ;  ’/,  Correlation  of  the  input  vector. 

Rx=Ca’*R=t=Ca; 

px=Ca'*R*Wq; 

wa( :  ,si+l)=wa( :  ,si)+mu*(px-Rx*wa(:  ,si)) ;  ’/,  Unconstrained  component  in  GSC. 

W(  :  ,si)=Wq-Ca*wa(  :  ,si)  ;  '/,  The  LMS  weight  estimate. 
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end 


7,  The  TDL  anteima  input  signal  vector  and  optimum  weights  value. 
si=0 ; 

for  i=Nl:N 
si=si+l 
for  n=l:K 
for  j=l:M 

ud(j+(n-l)*M,si)=inputd(i-Td*Rs*(j-l)-(n-l)) ; 
ui(j+(n-l)*M,si)=winputi(i-Ti*Rs*(j-l)-(n-l)) ; 
end 
end 
end 

u=ud+ui;  %  The  input  signal  vector. 

R=u*uV(N-N1)  7,  Correlation  of  the  input  vector. 

wa=inv(Ca’*R*Ca)*Ca’*R*Wq;  7.  Unconstrained  component  weights  in  GSC. 
Wopt=Wq-Ca*wa;  7.  Optimum  weight  of  the  LCMV-GSC  beamformer. 

7.  The  new  model  for  constraint  matrix,  input  vector,  and  LMS  weight  estimateW 
7o  in  the  frequency  hopped  environment. 
si=0 ; 
for  i=l:N 
for  j=l:K 
for  ii=l:M 

d(ii,l)=exp(-sqrt(-l)*2*pi*ffl(i)/2/(fl)*sin(thd)*(ii-l)) ;  7.  Array  response 

end  7o  vector. 
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C( : , j)=[zeros(M*(j-l) ,1) :  d;  zeros(M*(K-j) ,1)] ; 


'/,  The  constraint 


end  ‘/o  matrix 

f=[l  zeros(l,K-l)]  ’ ;  '/,  The  response  of  the  beamformer. 

Wq=C*inv(C’*C)*g;  */,  Wq  is  the  quiescent  weight  vector. 

[Ca]=block(C)  ;  */,  Blocking  Matrix 

for  n=l:K 
for  j=l:M 

ud(j+(n-l)*M,i)=inputd(i)*exp(-sqrt(-l)*2*pi*ffl(i)*(j-l)/2/(fl)\\ 
*sin(thd))*exp(-sqrt(-l)*2*pi*To*(ff (i))*(n-l)) ; 
ui(j+(n-l)*M,i)=winputi(i)*exp(-sqrt(-l)*2*pi*ff2(i)*(j-l)/2/(f 1)\\ 
*sin(thi))*exp(-sqrt(-l)*2*pi*To*(ff2(i))*(n-l)) ; 

end 

end 

u( : ,si)=(ud( : ,si)+ui( : ,si)) ;  %  The  input  signal  vector. 

R=u( :  ,si)*u( :  ,si)  ’ ;  ’/,  Correlation  of  the  input  vector. 

Rx=Ca’*R*Ca; 

px=Ca’*R*Wq; 

wa( : , si+l)=wa( : , si)+mu*(px-Rx*wa( :  ,si))  ;  Unconstrained  component  in  GSC. 
W(  :  ,si)=Wq-Ca*wa(  :  ,si)  :  "/.  The  LMS  weight  estimate, 

end 

B.3  Spatial  Smoothing  Technique 

’/,  Add  the  decorrelation  factor  in  the  LCMV-GSC  beamformer. 

R=zeros(M*K)  ;  7,  Initial  value  of  correlation. 

si=0; 

for  i=Nl:M 
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K1  =  5; 


'/,  The  smoothing  steps 


for  MN=1:K1 
for  j=l:K 
ti=0; 

for  ii=MN:M+Mn-l 
ti=ti+l; 

d(ti,  l)=exp(-sqrt(-l)*2*pi*f  i/2/(fl)=i'sin(thd)*(ii-l)) ; 
end 

C( : ,  j)=[zeros(M*(j-l)  ,  1)  ;  d;  zeros(M*(K-j)  ,  1)]  ;  ‘/.The  constraint 

end  '/,  matrix 

f=[l  zerosCl ,K-1)]  ’ ;  ‘/.  The  response  of  the  heamformer. 

Wq=C*inv(C’*C)*g;  7.  Wq  is  the  quiescent  weight  vector. 

[Ca]=block(C) ;  '/,  Blocking  Matrix 

for  n=l:K 
tt=0; 

for  j=MN:M+MN-l 
tt=tt+l ; 

ud(tt+(n-l)*M,si)=inputd(i-Td*Rs*(j-l)-(n-l)) ; 
ui(tt+(n-l)*M,si)=winputi(i-Ti=t'Rs*( j-l)-(n-l)) ; 
end 
end 

u(  :  ,si)  =  (ud(  :  ,si)+ui(  :  ,si))  ;  ’/.  The  input  signal  vector. 

R=R+u( : ,si)*u( : ,si) ’ ;  %  Correlation  of  the  input  vector, 

end 
R=R/K; 

Rx=Ca’*R*Ca; 
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px=Ca’=t=R*Wq; 


wa(  :  ,si+l)=wa(  :  ,si)+mu*(px-Rx*wa(;  ,si))  ;  '/.  Unconstrained  component  in  GSC. 

W(:  ,si)=Wq-Ca*wa(:  ,si);  */,  The  LMS  weight  estimate, 

end 


B.4  Adaptive  Filter 


inputd=[zeros(l,M-l)  input d] ; 
inputi=[zeros(l,M-l)  input i] ; 
input= [zeros input]; 
w=zeros(M,N+l) ; 
y=zeros(l ,N) ; 
e=zeros(l ,N) ; 
for  n=l:N-M 
for  j=l:M 

u(j ,n)=input(n+M-j) ; 

end 

y(n)  =  w(  :  ,n)  ’=t=u(:  ,n) ; 
e(n)  =  inputd(n+M-l)-y(n) ; 
w(:,n+l)  =  w(:,n)  +  mu*u( : ,n)*e(n) ’ ; 
end 


’/,  Put  a  leading  0  in  x 


’/,  Clear  tap  weights 
’/,  Observed  output 
■/,  Output  error 


7,  Input  vector 

7.  Output  point 
7.  Error 

7.  Filter  weights 
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